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ABSTRACT 

We present stellar population age and metallicity trends for a sample of 59 SO galaxies 
based on optical SDSS and NIR J & H photometry. When combined with optical g 
and r passband imaging data from the SDSS archive and stellar population models, 
we obtain radial age and metallicity trends out to at least 5 effective radii for most 
of the galaxies in our sample. The sample covers a range in stellar mass and light 
concentration. We find an average central light-weighted age of ~ 4 Gyr and central 
metallicity [Z/H] ~ 0.2 dex. Almost all galaxies show a negative metallicity gradient 
from the center out, with an average value of A[Z/H]/Alog(r/i? e ) = -0.6. An age 
increase, decrease, and minimal change with radius is observed for 58%, 19%, and 23%, 
respectively, for a mean age gradient of Aage/Alog(r/R e ) = 2.3 Gyr dex -1 . For 14 out 
of 59 galaxies, the light-weighted age of the outer region is greater than 10 Gyr. We 
find that galaxies with both lower mass and lower concentration have younger light- 
weighted ages and lower light-weighted metallicities. This mass-metallicity relation 
extends into the outer regions of our SO galaxies. Our results are consistent with the 
formation of SO galaxies through the transformation of spiral galaxy disks. Determining 
the structural component that makes up the outer region of galaxies with old outksirts 
is a necessary step to understand the formation history of SO galaxies. 

Key words: galaxies: elliptical and lenticular, cD; galaxies: evolution; galaxies: for- 
mation; galaxies: photometry 
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1 INTRODUCTION 

The unique morphological status of SO galaxies between 
gas-poor and gas-rich galaxies makes their origin a key 
ingredient of galaxy formation and evolution models. Nu- 
merous hypotheses suggest that at least a class of SO 
galaxies are by-products of secondary events in galaxy 
formation that transform spiral galaxies into SO galax- 
ies, rather than products of initial conditions <lickeMl985l : 
Bvrd fc ValtoneiJl99d: iMoore et aflll99rj ; lOuilis et aU200ol : 



Boselli et alj|2008h . These theories involve the removal of 



gas and subsequent halting of star formation and can in- 
clude the following formatio n mechan i sms: minor in terac- 
tions wit h other galaxies dlckel 1 19851 : iBarnesI (2002), ha- 
rassment (Mo ore et alJ|l99a). ram pressure and/or viscous 
stripping (Quilis ct al. 2000; Kronbcrgcr et al. 2008), stran- 
gulati on (iKawata fc Mulchaeyl [20081 ') . and cluster tidal ef- 
fects (|Bvrd fc Valtonenl ll99of ). Observational evidence for 
a higher fraction of SO-to-spiral galaxies in local clus- 
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ters as compared to higher redshift clusters certainly sup- 
ports the idea that SO galaxie s are primarily created 
through evolutionary processes l|Butcher fc Oemlerl 1 19781 : 
iDressler fc Gumj|l983h . However, obstacles to theories of 
a simple transformation from a spiral to a lenticular galaxy 
exist. For example, the bulge-to-disk ratios of SO galaxies 
appear to be on average larger than those of spirals in all 
density regimes (Dressier 1980J), thus creating problems for 
a straightforward disk-quenching scenario, although not for 
other scenarios such as merging. The structural dichotomy 
between massive and low mass SOs suggests diff erent forma- 
tion scenarios for the two mass groups (|van den Bergbiri994l : 
iBedregal et al.| [2008V In this paper we study the stellar pop- 
ulations in SO galaxies to better constrain their possible for- 
mation processes. 

Observations of radial stellar population (hereafter, SP) 
gradients in SO galaxies may help to distinguish between pos- 
sible scenarios for SO galaxy formation, since the dynamical 
events that lead to the production of SOs may produce or 
erase imprints on the galaxy's star formation history. For 
example, hierarchical models of galaxy formation predict an 
inside-out accretion of disk gas from the hot gaseous halo, 
resulting in ages that decrease outward in the galaxy (Fall & 
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Efstathiou 1980, but see Renzini 1995). On the other hand, 
formation mechanisms, such as those outlined above for the 
transformation of a spiral galaxy into an SO, could alter this 
tren d of age with radius through either recent star forma- 
tion llBvrd fc Valtonenlll990l;|Moore et al.ll 19961 ; "~ 
200a) or gas removal 



200C 


: 


Kronbcrgcr ct al 


2000 


: 


Kronberger et al. 



Quilis ct al 



Quilis et al 



2008). The latter process can alter 



age trends with radius because gas is preferentially removed 
from the outer disk, while for the former, this occurs through 
gas inflow. 

SP trends in the inner regions of SO galaxies have been 
widely studied. Spectroscopic studies suggest that the nu- 
clear regions of SO galaxi es are on average younger and more 



metal rich than bulges (IFisher et al.11199 1 



ISerra et al.ll2008l ); however. [Serra et all (120081 ) found this to 



: younger and more 
1 ISirchenkdl2006l ; 



be true only for Hi-poor SO ga l axies. Although the photo- 
metric studies of iPeletier et al.l (1 19991 ) found redder colors 
in the nuclear regions of SOs compared to further out, they 
attributed this to dust effects, as opposed to age differences. 
Most models for the transformation of spiral into SO galax- 
ies predict a resulting central gas concentration that leads 
to a young SP in the center of the SO rem nant (jlckel 1 19851 : 
iMoore et al"1ll996l ; iKronberger et~ai1l2008l ), consistent with 
spectroscopic studies (as cited above). Deri vations of the 
mean central ages yield an uncertain picture. IPeletier et alj 
i| 19991 ) examined the bulges of SO galaxies using Hubble 
Space Telescope images and found an a v erage bulge light- 
weighted age of 9 Gyr, while [Sil 'chenkol (2006) found ages 
between 4.8 and 8.3 Gyr. Both studies find high central 
metallicities. 

Even less is known about the age trends at larger radii 
in SO galaxies. Various authors have propose d that the disks 
of SO galaxies are younger than the bulg e |Caldwelllll983l : 
iBothun fc Gregdll990l : iNorris et all 12009). Others suggest 



comparable ages for the bu lge and disk (|Peletier fc Balceflsl 
Il999 : iMehlert et all 120031 ) while some report olde r ages 



in the disk compared to the bulge (Fisher et ahl 1 19961 : 



iTikhonov et al. | l2003l : iMacArthur ct al 



20041 ). The large 



variation in observed age trends may be due to differences 
in radial coverage, bandpasses studied, and the adopted SP 
models, if any, by the authors. Alternatively, it is possible 
that the disagreement between studies may be due to the 
heterogeneous nature of SO galaxies. Large samples with ac- 
curate photometry extending beyond a few effective radii 
are required to firm up SP and evolutionary trends in SO 
galaxies. 

In addition to establishing radial SP gradients in SO 
galaxies, one must assess whether the SP gradients and 
mean values are linked to global properties such as the 
galaxy's stellar mass, central light concentration, and envi- 
ronment. SO galaxies are proportionately more populous in 
the dense envi ronments of ri ch clusters than in low-density 
environments l|Dresslerlll980l ). suggesting that environment 
plays a role in SO formation and evolution. The mass of 
a galaxy is closely tied to its formation history. Hierarchi- 
cal assembly predicts that larger galaxy dark matter ha- 
los ar e formed through m ergers later in time than smaller 
halos l|Navarro et al.lfl995l ). However, semi-analytic models 
can produce "anti-hierarchical" star formation histories in 
a A Cold Dark Matter (ACDM) u niverse despite the h i- 
erarchic al assembly of t hese g alaxies (|De Lucia et al,| [2006). 
As well, iNeistein et al.l (|2006h showed analytically that anti- 



hierarchical evolution can be a natural outcome of bottom- 
up clustering. Such anti-hierarchical evolution is more con- 
sistent with actual observational trends of stellar popula- 
tions with galaxy mass for all Hubble types, including SOs 
l|MacArthur et al.l liooj . hereafter, M04). For example, an 
increasing fraction of SO galaxies with recent star formation 
has been observe d for decreasing luminosit i es (an indicator 
of stellar mass) l|Balcells fc Peletierl fl99i iPoggiant i et al.l 
l200lf ) . Additionally, a correlation of older stellar populations 
in SO galaxies with higher velocity dis persion (an indicator 
of dynamic mas s ) has been observed (jMehlert et al.l [20031 : 
ISil'chenkol [20061 : iBedregal et al.l 120081 ) . These studies sup- 
port different star formation histories for high and low mass 
SO galaxies, but have largely focused on the inner stellar 
populations of galaxies. 

The degree to which light is concentrated towards the 
center of the galaxy is also linked to a galaxy's formation 
history. Secondary galaxy evolution processes may increase 
the mass of the central regi on either through secular build - 
up or accretion of satellites (|Kormendv fc Kennicutt]|2004l ). 
On the other hand, the relative mass of the central com- 
ponent, or bulge, may influence the outcome of galaxy for- 
mation processes, such as those transforming spiral galax- 
ies. For example, harassment preferentially selects bulge- 
dominated galaxies to be stripped of their spiral structure 
and transformed into SOs, since disk-dominated galaxies will 
inste ad be shredded down to a dwarf system (|Moore et alj 
1996). Either way, the light concentration of the galaxy is 
closely tied to galaxy evolution processes and should thus 
correlat e with SP parameters. For instance. iMcDonald et ail 
(2009b) found a lack of low surface brightness galaxies with 
high concentration as well as a lack of high surface bright- 
ness galaxies with low concentration. Therefore, they sug- 
gest that galaxies with high concentration can only form 
through major mergers. Other studies relating the size of the 
bulge and the bulge-to-disk ratio (which is related to con- 
centration) to the evolution of SO galaxies find that galaxies 
with la rger bulges are redder a nd have smaller color gra- 
dients dBalcells fc Peletierl II994T) . How e ver, a t intermediate 
redshift (0.73 < z < 1.04). IKoo et al.1 l|2005l ) found no dif- 
ference in bulge colors between disk and bulge-dominated 
galaxies, suggesting that B/D ratio is poorly correlated with 
stellar populations. 

Despite recent progress, studies of radial SP trends in 
SO galaxies as a function of the global properties of mass 
and concentration are still open-ended and have seldom ex- 
plored the regions beyond 1 — 2 R e . Most efforts have indeed 
focused on age and chemical composition within one effec- 
tive radius. In this paper, we use deep optical/near-infrared 
(NIR) color gradients to explore age and metallicity gradi- 
ents in SO galaxies out to typically 5 R e , and with a large 
enough sample to cover a substantial range in galaxy mass 
and light concentration. The focus of this paper is largely 
on empirical radial trends in stellar populations; separate 
bulge and disk trends that result from model decomposi- 
tions of the galaxy image will be presented in a forthcoming 
publication (Chamberlain et al. 2011, in prep). Our sample 
of 59 SO galaxies is presented in jj2] while our observations 
and data reduction techniques are discussed in ij3] and ^4] 
respectively. Measurements of global galaxy properties are 
presented in ij5]and the SP models are described in i|6] Our 
results for the age and chemical composition trends in SO 
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galaxies are presented in 33 The implications of our results 
are discussed in SJH an d a summary is presented in iJU A 
distance of 16.5 Mpc is used for all Virgo cluster galaxies 
(Mei et al. 2007). 



2 SAMPLE 

We selected from the UGC catalog (Nilson 1973) all galax- 
ies with SO, SOa and SOB morphologies and blue Galactic 
extinction < 0.5 yielding a catalog of 1088 objects. This 
sample was further restricted to areas of the sky covered by 
the Sloan Digital Sky Survey l|York et all I2000L hereafter, 
SDSS) bringing the sample to 542 galaxies. Finite telescope 
time limited us to NIR data for only ~ 15% of the SDSS 
sub-sample; the sample was selected to cover a represen- 
tative range of mass and light concentrations. Eliminating 
edge-on galaxies and galaxies with significant spiral struc- 
ture or overlapping foreground stars further restricted our 
sample to 59 objects as described below. While this UGC- 
derived sample includes some range in environment (see 
SJSJ), the higher density environments are missed. Thus, for 
greater coverage, we include observations of eight low-mass 
Virgo cluster galaxies from the sample of iMcDonald et aU 
(2010) (hereafter M10); these were observed with the same 
equipment as our sample (see £j3]). However, these eight 
galaxies all have low mass and low concentration; see jJH 
for more details. In the electronic version, we show SDSS 
color images (using g, r, and i bands) for our full sam- 
ple. These images were obtained from the SDSS website: 
|http://cas.sdss.org/dr7/en/tools/chart/navi.asp| 

Our sample excludes galaxies with spiral structure and 
highly disturbed structure, as gauged by visual inspection 
of the SDSS images. Galaxies with tidal tails or faint spiral 
structure were retained, but classified as "transition" galax- 
ies due to the emergence or fading of their spiral structure. 
Examples are shown in the color images available in the 
electronic version, and are noted by a "T" above the image. 
Twelve of our 59 sample galaxies are labeled as "transition" . 
We will see in £17.31 that the putative structural transition 
correlates with a transition in SP properties. We note that 
some SO galaxies left in our "featureless disk" sample are 
likely to reveal spiral structure when examined with higher 
resolution images and our sample may also contain elliptical 
galaxies misclassified as SO. 

Highly inclined galaxies with i > 75 deg, and galaxies 
that fall on the edge of the SDSS field of view, were also 
pruned from the sample in order to measure well-defined 
surface brightness (SB) profiles. 

Our final sample is neither statistically complete nor 
randomly selected, but tailored to cover a broad range in 
global properties such as mass, concentration, and local en- 
vironmental density (see The galaxy sample is presented 
in Table 1. 



3 OBSERVATIONS 

Our SP analysis relies on optical and NIR photometry of 59 
SO galaxies. The NIR observations (J- &/or H-band) were 
collected with the ULBCam at the University of Hawaii's 



2.2-m telescope and the archival optical imaging (g & r) is 
from SDSS. 

3.1 Optical Observations 

Optical g & r images with a scale of 0.396 arcsec/pixel were 
extracted from the SDSS/DR5 archives. The u and z bands 
were left out considering their low signal-to-noise. We did 
not use the i-band since it does not improve the discrim- 
inating power for age and metallicity in SP models, when 
combined with g, r, and J or H bands. 

3.2 Near-IR Observations 

Near-IR images for the 59 SO galaxies in our sample were 
obtained with the ULBCam at the University of Hawaii's 
2.2-meter telescope on Mauna Kea in April 2005, 2006 and 
2007 and March 2008. A total of 29 galaxies in our sample 
were observed only in the H band, 19 were observed only 
in the J band, and 11 galaxies were observed in both the J 
and H bands. The ULBCam image scale is 0.25 arcsecond 
per pixel. A maximum single exposure time of 40 seconds 
was used to maximize the sky flux whilst keeping within the 
detector's linear regime. The total integration time depends 
on SB and ranges from 8 to 24 minutes. A standard dither- 
ing script minimized the resampling of bad pixels. Only the 
cleanest of the four 2048x2048 arrays was used, resulting 
in an 8.5' by 8.5' field of view, which was large enough to 
properly assess the sky background levels. Our NIR flux cali- 
bration makes use of reference stars in the target galaxy field 
from the Two Micron All-Sky Survey (|Skrutskie et alJbood . 
hereafter 2MASS). Thus, no standard star observations for 
photometric calibration were necessary (see i|4.1[) . A more 
detailed description of ULBCam data taking procedures is 
provided in M10. A log of observations is presented in Table 
2. 

In the remainder of the paper, we refer to the galaxies 
observed in the J-band as the "J-band sample" and those 
observed in the H-band as the "H-band sample". Eleven 
overlapping galaxies are included in both samples. 



4 DATA REDUCTION 
4.1 Basic Reductions 

SDSS images are already processed for fiat-fielding, bias sub- 
traction, and cosmic ray rejection. The photometric calibra- 
tion of the light profiles ( t|4.6|) used the photometric zero- 
points provided in the SDSS/DR5 library. 

Basic reductions of the UH NIR data, which include 
flat-fielding, stacking, bad pixel rejection, geometric distor- 
tion corrections, and flux calibration, were applied to the 
ULBCam data using the XVISTA software packagtQ- The 
data reduction procedures follow the prescription of M10. 
The NIR flux calibration of light profiles used infrared stel- 
lar photometry from 2MASS for stars in the ULBCam tar- 
get galaxy fields. A mean offset between our brightnesses 
and the 2MASS J or H-band brightnesses was calculated 

1 http://astronomy.nmsu.edu/holtz/xvista/index.html 
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and used for the zero-point calibration. This method of flux 
calibration enabled us to calibrate the photometry of our 
science images at each pixel, independent of airmass vari- 
ations and transparency conditions. Additional information 
on flux calibration methods and background stability can be 
found in M10. The photometric zero-point and the standard 
deviation about the mean in the calibrating stars are noted 
in Table 2. If fewer than three 2MASS stars were found in a 
target field (indicated by a zero in the last column of Table 
2), the calibration used is the average of the previous and 
next exposures. The quoted standard deviation for these ob- 
servations is a standard deviation of all calibration values, 
that is ± 0.056 mag. 

SBs have been corrected for Galactic foreground extinc- 
tion in each bandpa ss using the reddening values, A\, of 
ISchlegel et all (|l998l ) and assuming an Rv = 3.1 extinction 
curve. 

4.2 PSF Matching 

Wavelength-dependent image blurring by the atmosphere af- 
fects inner galaxy colors. Thus, we measured the 2D seeing 
point spread function (PSF) to correct galaxy color profiles. 
The PSF full width at half maximum (FWHM) is measured 
from the final stacked image in each bandpass. The stars 
that are identified for PSF measurements in the g & r bands 
are individually sky-subtracted and fit with a 2-dimensional 
Gaussian function from which FWHMs are derived. On aver- 
age, 50 stars are used for PSF measurements in each stacked 
SDSS image. The final PSF per image is the median of all in- 
dividual star PSFs. For NIR images, the final PSF per image 
is the average of 5 individual PSFs. Images from each band- 
pass are degraded to the PSF of the worst seeing bandpass 
(typically the g image) by convolving with a 2D gaussian of 
appropriate FWHM. For additional precaution, we apply an 
inner radial cutoff of roughly 2 seeing disks (3 arcseconds) 
to the color profiles. The central region in all other analy- 
ses of colors is averaged within at-least the inner 2 seeing 
disks. PSF measurements have also enabled the identifica- 
tion of foreground stars that are then removed from the 
galaxy light. 

4.3 Sky Measurement 

Careful sky subtraction is crucial for the accurate analysis 
of deep SB profiles. Sky subtraction follows slightly different 
approaches for SDSS and NIR images. 

We have used the sky value provided with each SDSS 
imag60. For the images where the sky value was not available 
(~ half our sample) , we used the lowest of the sky measure- 
ments measured by us from either the entire SDSS image 
or from the four image corners. The latter scheme gives sky 
estimates within 0.5% of the SDSS values (M10). 

For the NIR data, the sky is measured in four rectangu- 
lar boxes along the perimeter of the field of view and then 
averaged together to give the mean sky value for that image. 
Each of the four sky boxes has a typical size of ~ 75 x 500 
arcseconds. The typical deviation of sky levels amongst the 
four boxes for the NIR data are .004% of the sky value. 

2 http:/ /www. sdss.org/dr5/algorithms/sky.html 



4.4 Surface Brightness Profile Extraction and 
Error Estimates 

SB profiles were extracted by fitting elliptical isophotes to 
the r-band galaxy images. The XVISTA command, PRO- 
FILE, enables this operation based on a generalized non- 
linear least-squares fitting routine. For these fits, ellipticity 
and position angle are allowed to vary but the galaxy center 
is fixed. The isophotal solutions based on the r-band images 
were then applied uniformly to the images in the g, J, and 
H band images. This ensures that color gradients are com- 
puted from the same matching isophotes. Further details 
about isophotal fitting and profile extraction, including de- 
tails on p r ofile d epth and signal-to-noise ratios, are given in 
ICourteaul (|l996l ) and M10. 

The effect of a systematic sky error on galaxy photom- 
etry naturally increases with galactocentric radius. To esti- 
mate this effect, we have recalculated the SB profiles and 
color profiles using sky values adjusted to sky = sky or i g 
±la sya , where a ays is the systematic sky uncertainty. This 
is calculated as the standard error in the mean of the sky val- 
ues (as described in ij4.3[l from four boxes along the perime- 
ter of the image. We determine the sky effect at each radius 
along the SB profile as half the difference in the SB calcu- 
lated with the sky value set at sky + = sky or i g + a sys and 
sky~ = skyorig - o- 3ys . Our final computation of the photo- 
metric error at each radius is the sum in quadrature of the 
statistical error from profile fitting (this naturally includes 
statistical errors in the sky background) and the systematic 
sky error estimate. 

4.5 Outer Radial Cutoff 

The low NIR galaxy SB, relative to the bright NIR sky 
background, make an accurate determination of the sky 
background crucial for extracting reliable surface photom- 
etry. This is demonstrated in Fig. [1] for a sample galaxy 
(UGC4737). We define the outer radial cutoff as the point 
where the sky error envelopes deviate from each other by 
more than 0.4 mag arcsecond~ 2 , which is usually where the 
statistical error per point is comparable to the typical point- 
to-point variations. Our NIR profiles are intrinsically shal- 
lower than SDSS profiles and the common outer truncation 
radius was thus determined using J and H profiles for all but 
UGC 10158 where the optical bandpasses are more sensitive 
to sky errors. In the few cases where the statistical error 
from profile fitting (for any of the optical or NIR bands) 
consistently exceeds 0.1 mag arcsecond -2 at large radii, the 
radius at that SB threshold point would be used as the outer 
cutoff instead. An example of the outer cutoff, based on an 
H-band profile, is shown in Fig. [T]as the dashed vertical line. 
The radial cutoff used for the g and r bands matches the 
NIR band that is being used for analysis. The outer radial 
cutoffs for each galaxy in the J and H-bands are given in 
Table 3. 

4.6 Surface Brightness and Color Profiles 

Optical (g & r) and NIR (J and/or H) SB profiles are ex- 
tracted for all of our galaxies following the prescription out- 
lined in t]4.4l Fig. [2] shows an example surface brightness 
profile for UGC 4869. The entire collection of SB profiles is 
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shown in Fig. [3] The g, r, J and H bandpasses are repre- 
sented by blue, magenta, red and black data points. The J 
and H profiles are displayed for the 11 galaxies observed in 
both passbands. Color profiles for g-r and r-H versus radius 
are shown in Fig. [4] In figures where both J and H profiles 
are displayed, the cutoff for the g and r passbands is deter- 
mined by the larger of the J and H-band cutoffs. The SB and 
color error bars at each radius reflect the ± la statistical er- 
rors only. These errors tend to be small, generally smaller 
than the point size, except at the largest radii. Interior to 
the outer radial cutoff, SBs and colors are not plotted for 
radii where the statistical error exceeds 0.1 mag arcsec -2 , 
which is equivalent to setting very low weight for these pix- 
els. Gaps in SB and color profiles, such as that seen for 
UGC04737 around 25", are due to overlapping foreground 
stars or galaxies. Since these foreground objects have been 
masked during ellipse fitting, which is identical to setting 
very low weight for these pixels; they do not affect the pro- 
file shape. 

We compute effective radii and total magnitudes in the 
g, r, J, and H bands. To determine the total magnitude, we 
have extrapolated the SB profile outwar d by fitting the o uter 
profile with an exponential function. See lCourteaul (|l996T ) for 
more details. The effective radius, R e , is the radius that con- 
tains 50% of the total extrapolated light. F or the 5 galaxie s 
with a plateau in their SB profiles (Type II; iFreemanll 197oT ) . 
no extrapolation is performed. Total magnitudes, mean col- 
ors and Re for each galaxy in our sample are given in Table 
4. 

The lower axis of the SB and color profiles in Figs. [3]- 
[4] shows the radius in terms of the r-band R e of the galaxy. 
Most profiles extend past 5 R e . It is also apparent that our 
sample galaxies cover a range in profile shapes, suggesting 
a variety of contributions from galactic components. Dips 
and plateaus in the SB profiles, such as those seen for UGC 
4596, are due to structure in the galaxy, usually rings, bars, 
or spiral structure, and/or extinction by dust. 

4.7 Radial Binning 

To increase the signal per color bin, we average colors in 
six radial bins, scaled by the r-band half light radius. The 
binning scheme is noted in Table 5. If the statistical SB error 
at any radial point exceeds 0.1 mag arcsec -2 , that radial 
point or region is excluded from the analysis. We also ensure 
that each binning region contains at least 5 radial points, 
otherwise that region is excluded from further analysis. 

Each galaxy is also subdivided into an inner and outer 
radial region. The separation between the two regions is cho- 
sen at the most prominent inflection point in the SB profile. 
We find a well-defined change in the slope of the SB pro- 
file for ~ 50% of the sample, with the majority having a 
break between 0.8 and 1.2 R e . We thus label inner regions 
as those within 0.8 R e and outer regions as those beyond 
1.2 R e for all galaxies in our sample. Fig. [2] shows an ex- 
ample of the inner and outer regions with dashed lines at 
0.8 and 1.2 R e . Note that our distinction of inner and outer 
regions is independent of model fitting of galactic compo- 
nents (i.e. disk and bulge). Analysis of stellar populations 
of discrete galactic components will be discussed elsewhere 
(Chamberlain et al. 2010, in preparation). Here we focus on 
a model-independent analysis of the SB profiles. For clarity, 



when discussing the central-most colors, we refer to 'central' 
as the binning region with r < 0.5 R e and 'inner' for r < 0.8 

Re- 



5 DETERMINATION AND RANGE OF 
GALAXY PROPERTIES 

Our analysis of SP gradients relies on correlations with 
global galaxy properties, such as stellar mass, light con- 
centration and environment. Prescriptions for determining 
these properties are given below and the values for each 
galaxy in our sample are listed in Table 6. 

Total stellar masses for our galaxies were calculated 
from total g and r magnitudes (described in £1-4.61) a nd us - 
ing the mass-to-light ratio prescriptions of lBell et al.l l|2003h , 
unless otherwise notecjf]. Distance estimates, which are re- 
quired to calculate the mass and physical size of the galax- 
ies, were corrected for Virgo flow and the Great Attractor 
as provided by the NEE0 for all galaxies other than Virgo 
cluster galaxies. The median stellar mass of our sample is 
1 x 10 11 Mq. Fig. \S\ shows a slight bias towards high mass 
SOs. 

Model-dependent bulge-to-disk ratio estimates may 
carry large systematic errors due to the sub j ectiv- 
ity of profile fitting functions |MacArthur et al] 120031 ; 
iMcDonald et al]|2009al ). Alternatively, the galaxy light con- 
centration parameter g ives a non- parametric indication of 
the bulge-to-disk ratio (Kent 1985). It is computed as: 

C 2 8 = 51og(r 8 o/r 2 o) 

where rso and r 2 o are the radii that enclose 80% and 20% 
of the total light, respectively. The total magnitudes and 
subsequent 80% and 20% radii have been calculated as de- 
scribed in i]4.6l Concentration values for our sample galaxies 
range from C28 = 2.8 to 5.5 with a median value of 4.7 as 
shown in Fig. For reference, a pure exponential disk has 
C28 ~ 2.8. 

Our definition of local environment uses a three- 
dimensional number density based on the mean distance of 
the six nearest neighbors. To constr uct the density fie ld, we 
use the Updated Zwicky Catalog l|Falco et al ] Il999l , here- 
after UZC) which is 95% complete to a limiting magni- 
tude of mz, = 15.5 mag. A three-dimensional position in 
a Cartesian coordinate system of each galaxy in our sample 
is then assigned, based on its sky position and recessional 
velocity (using a value for the Hubble constant of H = 75 
km s _1 Mpc _1 ). The local number density of our SO galax- 
ies is calculated by taking the number of objects contained 
within a sphere of radius equal to the mean distance of the 
six nearest neighbors in the density field and then dividing 
by the physical volume of the sphere. We note several sources 
of uncertainty in our environmental density estimate. First, 
the magnitude limited catalog samples the galaxy luminosity 
function to different levels at different heliocentric distances. 

3 We have also calculated masses using iPortinari et al. 1 mil) 
and find no difference in our results 

4 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration 
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To correct for this effect, our densities are multiplied by a 
luminosity function (LF) correction factor, which is the ratio 
of the integrated observable LF (i.e., the number of galaxies 
per Mpc _3 brighter than the limiting magnitude of the sam- 
ple) at 3000 km s _1 to the integ rated LF at the galaxy's 
redshift. We use a Schechter LF (|Schechterl ll976T) derived 
from the UZC, with a=-1.0 and M*=-18.8 (|Marzke et"afl 



1994). However, for objects with cz > 9000 km s 
rection factors become unreliable (greater than factors of 
3). Fortunately, only 8 out of our 59 galaxies fall beyond 
this redshift range. We also consider that galaxies with cz 
< 3000 km s" 1 (17 galaxies in our sample) have large pe- 
culiar velocities that yield uncertain line-of-sight distances. 
Finally, in clusters, the high internal velocity dispersion will 
bias number densities towards lower values. The distribution 
in environmental densities for our sample covers the range 
from -2.2 log Mpc" 3 to 1.23 log Mpc" 3 with a median 
value of -0.67 log Mpc -3 . Our sample is underrepresented 
in the cluster regime; our 8 Virgo cluster galaxies sample 
the densest environment. The environmental density calcu- 
lations were kindly provided by Jesse Miner. 

Fig. [6] demonstrates the broad parameter space covered 
by our sample. Black crosses and green dots refer to galaxies 
in the lowest and highest density environments respectively. 
Although our SO galaxies have mostly low concentration for 
low mass galaxies, SP trends in mass and concentration can 
still be identified. However, the combination of a statisti- 
cally incomplete sample regarding environment (e.g., all of 
the observed lowest mass galaxies are in the Virgo cluster 
high density environment) and a scarcity of high density 
environment galaxies makes it difficult to draw conclusions 
regarding the effect of environment on age or metallicity. 

Galaxies typically occupy two distinct regions in color- 
stellar mass s pace; the so-called " red sequence" and the 



"blue cloud" Jstrateva et alj l200ll; iKauffmann et al.l [20o3 : 
iBell et alJ 12004 iKannappan et all l2009eft . Although SO 
galaxies generally fall into the red sequence, studies have 
shown that thi s morphology- color corre spondence fails for 
low mass SOs ( Kannap pan et all 12009 bT) . Fig. [7] shows u- 
r color (Petrosian ma gnitudes from SPSS) versus stellar 
mass calculated using iPortinari et~ai] (|2004h . The dotted 
line shows the boundary between red an d blue sequences, 
as reported bv iKannappan et alJ l|2009bl ) but with a mass 
offset facto r of 1.8 to account fo r mass scale differences. 
We use the Portinari ct al. (2004) mass transformations to 
match with lKannappan et al.l (|2009bl V All but one SO galaxy 
(UGC09003) in our sample fall in the red sequence. 



6 STELLAR POPULATION MODELS 

A comparison of the observed galaxy colors with those pre- 
dicted from SP synthesis models covering a range in age and 
chemical composition allows for the determination of light- 
weighted mean ages and metallicities. The combination of 
a primarily age sensitive color (e.g., g-r) with a primarily 
metallicity sensitive color (e.g., r -H) provides a separation 
in age and metallicity. We use the lCharlot fc Bruzuall l|2010l. 
private communication, hereafter CB10) simple stellar pop- 
ulation (hereafter, SSP) model with a Salpeter initial m ass 
function and STELIB libraries |Le Borgne et all 120031 ') for 
our analysis. Since only a single age is attached to any SSP 



model whilst the observed light likely results from multiple 
coeval populations, the derived ages of the stellar popula- 
tion contributing the integrated light are best understood 
as a light-weighted mean age rather than a mass-weighted 
age. From here on, we drop the "mean" and refer to these 
as light-weighted ages and metallicities to avoid confusion 
with mean ages that are radially averaged. 

Ages and metallicities are presented in Fig. [8] for the 
UGC10391 color-color diagram. For each galaxy, the central 
bin, subsequent radial bin, and consecutive radial bins are 
denoted by a green star, small filled circle, and a solid blue 
line, respectively. Overplotted is a CB10 SSP model grid. 
Red dashed lines represent model lines of constant age in- 
creasing, left to right, from 0.8 Gyr to 13.8 Gyr. Blue dotted 
lines represent model lines of constant metallicity increasing, 
bottom to top, from [Z/H] = -2.2 to +0.5. The error bars 
for each radial bin represent the standard error in the mean 
based on the scatter in color of the radial points within the 
designated bin added in quadrature to the sky effect (see 
34T3J ■ This example shows a galaxy whose light-weighted 
metallicity and light-weighted age decrease and increase re- 
spectively from the center of the galaxy outward. The galaxy 
gets bluer in r-H at larger radii and keeps a roughly constant 
g-r, thus crossing over lines of constant age with increasing 
radius and indicating older ages in the outer regions of the 
galaxy. 



6.1 Star Formation Histories and Model 
Uncertaint ies 

Since our analysis of SP trends in SOs hinges on the reliabil- 
ity of derived ages and metallicities, we must understand the 
effects of SP model differences on our results. The SSP mod- 
els that we use for our analysis predict the evolution in col- 
ors (and spectrum) of a coeval population of stars with the 
same chemical composition and specified initial mass func- 
tion. While this is clearly an over-simplification of actual star 
formation histories in SO galaxies, it represents a straight- 
forward way to obtain a light-weighted age and metallicity. 
However, we must test how our analysis could be affected 
by the use of more complex star formation histories. Since 
an SSP is one extreme star formation history (equivalent to 
a single burst), we compute models for the other extreme, 
a constant star formation history with a quenching of star 
formation at various ages, and compute models for the inter- 
mediate case of exponentially declining star formation. We 
thus consider populations that are composed of a superposi- 
tion of SSPs, born at different epochs. Using the csp-galaxev 
prograrrjfl, we use an SSP model with a Salpeter initial mass 
function and fixed metallicity, convolved with the given star 
formation history (exponentially declining and constant). 
We compute the convolved models for a range of metallici- 
ties and time constants, r, (for the exponentially declining 
star formation models) or star formation truncation times 
(for the constant star formation history models). In both 
cases, we set the age of the galaxy, or time that star forma- 
tion began, to 13 Gyr. We show a color-color diagram based 



5 Part of the iBruzual fc Chariot! [1200311 software release: 
http: / / www2.iap.fr/users / charlot/bc2003/ 
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on the two sets of models in Fig. [9] The constant star for- 
mation history is shown in blue, the exponentially declining 
star formation history in red, and the SSP model in black. 
The same SSP model grid is used throughout the paper. 
The meaning of the "age" of an SP for each set of models 
is different. For the exponentially declining models, we plot 
lines of constant r. For the constant star formation history 
models, the age is represented by truncation times, or how 
long the star formation has lasted since 13 Gyr. The model 
grid edges are similar between the three model sets. For ex- 
ample, the constant star formation model with a truncation 
time of 0.1 Gyr and the exponentially declining model with 
a r of 0.1 Gyr are similar to the SSP model with an age 
of 13.8 Gyr, all lying and nearly overlapping near the right 
edge of the grids. The general shape of the lines do not vary 
with star formation history, except for the variation of the - 
0.3 metallicity line at young ages which is small compared to 
the large metallicity gradient that is observed in our sample. 
Thus, an analysis based on relative age trends with radius 
will be robust with respect to star formation histories. 

Uncertainties in SP synthesis modeling have received 
notable attention (Trager et a l. 2000; Schi avon et al.l [20021 ; 
IShaplev et ail 12003 ; IConrov et al.ll2009l ). Theoretical model 
uncertainties in the age and metallicity zeropoint can be 
caused by errors in the calibration of ages and metallici- 
ties from globular clusters. When two populations differ in 
metallicity as well as in age, and if the model zero point 
errors are metallicity dependent, differential ages between 
metal-poor and metal-rich populations are less secure. In- 
deed, there appears to be a strong metallicity gradient with 
radius for most of our galaxies. Therefore, our derived age 
trends with radius are subject to metallicity- dependent zero 
point issues and we take this into consideration when ex- 
amining our results. Another concern is the difference in 
elemental abundance ratios between the population being 
studied and the stars from which the models are based. 
Massive SO galaxies are known to have a higher abun- 
dance of alpha elements, relative to iron, as compared to 
solar neighborhood stars dNelan et al.ll2005l ; pTh omas et al.l 
l2005l ; ISerra et al.ll2008l ; lBedregal et al.ll2008t ) . Thus when us- 
ing stellar population models that do not account for non- 
solar abundance ratios, derived ages and metallicities are af- 
fected by errors in the theoretic al stellar evolutionary tracks 
l|Schiavonll20"07l ; iLee et all [20091) . Uncertainties in the input 
model parameters, such as the effective temperature of the 
isochrones, giant stars or binary stars in the luminosity func- 
tion, our understanding of late evolutionary phases, as well 
as the theory for convection and the effects of rotation and 
diffusion, may induce additional, unknown systematic er- 
rors. The treatment of advanced stages in evolution, such 
as the Thermal-Pulsating Asymptotic Giant Branch phase 
(TP-AGB), has received much attention lately. TP-AGB 
stars are extremely bright and dominate the NIR light of 
a galaxy following a burst of star formation, but represent a 
challenge for theoretical models due to the combined effects 
of thermal pul ses, changes from heavy element dredge- up, 
and mass loss l|Bruzual fc Charlotl [20031 ; lMarastonll2005l ). 

The significant range in predictions between different 
SSP models can lead to great differences in the interpre- 
tation of galaxy colors. Although it is difficult to estimate 
uncertainties in the derived SP parameters (like age and 
metallicity), the comparison of available models suggests 



how re liable they are. The CB10 model is a revised version 
of the iBruzual fc Chariot! (|2003l . hereafter, BC03) models, 
which includes a new prescription for TP -AGB evolution of 
low an d intermediate mass stars following iMarigo fc Girardil 
(2007), and uses tracks and isoc h rones with updated i nput 
physics fromlBertelli et all (|2008h .l Eminian etall <|2008l ) and 
iMacArthur et all (|201oF have highlighted s gnificant changes 
with the new Bruzual and Chariot models in NIR model 
colors for inter mediate populat ions. As well, the SP syn- 
thesis models of lMarastonl (|2005l ') different from the Bruzual 
and Chariot models by using a fuel-consumption-based code. 
These models place a particular emphasis on the TP-AGB 
phase and the effect on the model colors h as been demon- 
strated in the literature l|Tonini et al .1120091 ). 

We c ompare the origin al BC03 SSP models t o the SSP 
model s of lMarastonl (|2005l ) and CB10 in Fig. [TO] iMarastonl 
(2005) and CB10 SSP model grids are shown in the left and 
right panels, respectively. The BC03 model grid is shown in 
each panel (in black) for reference. All three models cover 
similar ages and metallicities (see figure caption for details), 
but the lowest age plotted here (0.8 Gyr ) is not available 
for the lowest metallicity (-2.25 dex) in the lMarastonl (120051) 



Maraston 
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model . The CB10 models are a closer match to the 
(2005) models than the earlier BC03 version. Both 
(2005) and CB10 models are redder in r-H, which pro- 
vides a closer match t o observations of star forming galaxies 
l|Eminian et al.l [20081 ). Despite large differences in color at 
low ages, we find in all models that the lines of constant age 
are tilted in a similar direction at all ages. We will return to 
this key point in H7.2I For younger ages ( ~ 3 Gyr), the lines 
of constant metallicity in the CB10 and IMarastonl (|2005l ) 
models change slope (r-H colors decrease along lines of con- 
stant metallicity for increasing age) and lines of constant 
age vary slightly in their dependence on r-H colors. We will 
discuss the model differences on colors in SJ7] 



6.2 Extraction of Ages, Metallicities, and 
Gradients 

For our analysis, we use the CB10 models with an SSP star 
formation history. We derive light- weighted ages and metal- 
licities for each radial bin and for the inner and outer ra- 
dial regions (as denned in >14.7|) by fitting CB10 SSP mod- 
els to the galaxy colors. We compute a finely-spaced CB10 
SSP model grid by interpolating linearly between the SSP 
metallicities and between the finely spaced ages provided by 
CB10. From the g-r and r-H (or r-J) colors, we determine 
the ages and metallicities from the model that minimizes the 
difference between model and observed colors, calculated in 
quadrature. Due to the uncertainties in extrapolating model 
colors to larger metallicities and ages and the convergence 
in colors at large ages, the age and metallicity of a galaxy 
are set to a maximum of 13.8 Gyr and +0.5 dex, respec- 
tively. While the CB10 models extend to 20 Gyr, we chose a 
maximum age of 13. 8 Gyr for consistency with the current 
age of the universe (|Hinshaw et al.ll2009l ) (the colors of 20 
Gyr and 13.8 Gyr models are similar enough not to affect 
our analysis). 

The uncertainties in the derived ages and metallicities 
are determined via a Monte Carlo method based on M04. 
Two-hundred realizations of the model fits were performed 
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for each radial bin and inner and outer regions with the col- 
ors for each realization drawn from a Gaussian distribution 
of the errors in each color. The standard error in the mean 
is based on the scatter in color of the radial points within 
the designated bin or region added in quadrature to the sky 
effect. Ages and metallicities for each bin correspond to the 
mean of the ages and metallicities computed from each re- 
alization. The positive and negative errors in the derived 
ages and metallicities for each bin are then taken as the in- 
terval containing 68% of the 200 Monte Carlo realizations. 
If more than 5% of the realizations produced either an age 
greater than 13.8 Gyr or a metallicity greater than +0.5, 
we set the error to zero, indicating the lack of a measurable 
error. We note that model fitting errors do not reflect un- 
certainties in the model itself (see N6.ip . Ages, metallicities, 
and their uncertainties are listed in Table 3 for the inner 
and outer regions of all galaxies in our sample. Tables of 
ages and metallicities at each radial bin are available upon 
request. 

There are 11 galaxies with available J and H-band im- 
ages, enabling a partial consistency check on our extracted 
ages and metallicities. In Fig. 1111 we show the derived J- 
band ages versus the derived H-band ages for both the in- 
ner and outer regions. A Kolmogorov-Smirnov two-sample 
test (hereafter, KS) on ages for the 11 galaxies in common 
reveals a probability of 81% that J and H-band derived ages 
are drawn from the same population for both inner and outer 
regions. 

Age and metallicity gradients are determined from a 
linear least-squares fit to the age and metallicity data in the 
second, third, and fourth radial bins. We avoid the central 
radial bin due to the greater likelihood of dust contami- 
nation and seeing blur; the fifth and sixth bins are avoided 
due to their greater errors from systematic sky uncertainties. 
The age and metallicity gradients for each galaxy are listed 
in the electronic version. Because the extracted ages and 
metallicities were set to a maximum of 13.8 Gyr and +0.5 
dex respectively, when an age (outermost bin) or metallicity 
(innermost bin) is set to one of those values, the gradient is 
considered a lower limit. 



7 RESULTS 

We wish to constrain the star formation and chemical enrich- 
ment histories of SO galaxies by mapping light- weighted ages 
and metallicities from the center out to large radii. We first 
analyze the central ages and metallicities. While sky back- 
ground removal is clearly the limiting obstacle to accurate 
photometric colors at large radii, it is a minimal problem in 
the central region. In addition, there exists numerous age de- 
terminations for the centers of early-type galaxies, including 
SOs, in the literature, thereby enabling a useful comparison 
of modeling procedures. On the other hand, central regions 
of early-type galaxies, which may include SOs, can be "con- 
taminated" by recent star formation episodes that may be 
non-r epresentative of the mean age of the rest of the galaxy 
fe.g.. Ide Jong fc Davieslll997l : ISanchez-Blazquez et al.|[2007l ; 
iMacArthur et alj|2009l ) Therefore the mean age derived for 
the central region may depend on the aperture size for the 
age determination. Moreover, PSF differences between dif- 
ferent passbands make the interpretation of surface photom- 



etry within the central regions of galaxies problematic. De- 
spite these complications, we first compare our central ages 
and metallicities with corresponding values from the litera- 
ture. An analysis of the galaxies' outer regions will follow. 
We note that our central regions all correspond to a radius 
of at least 3 arcseconds to alleviate PSF matching issues. 

7.1 Central Ages and Metallicities 

Central g-r colors versus r-H colors are plotted in Fig. ll2l for 
all galaxies in the H-band sample relative to a SSP model 
grid, where the central bins (r < 0.5 R e ) are designated by a 
green star. It can be seen that most of our galaxies have cen- 
tral colors that are concentrated in regions of intermediate 
age and high metallicity. In Fig. [13] we plot histograms of the 
extracted central (r < 0.5 R e ) ages and metallicities for the 
J- and H-band samples. The 11 galaxies with J- and H-band 
data yield similar age and metallicity histograms, therefore 
we shall merge the two data sets. However, we note non- 
uniformities in the J-band backgrounds that do not exist at 
H. Thus, when combining all the galaxies in our sample, the 
quoted ages and metallicities for galaxies that have H-band 
imaging will be based on that band alone. 

The distribution of light-weighted central ages for the 
combined J- and H-band data is fairly symmetrical with 
a mean age of 3.6 ±1.8 Gyr (the error is an rms scat- 
ter). The median central age is 3.3 Gyr. Thus, fairly re- 
cent star formation episodes have occurred in the cen- 
tral regions of a majority of SO galaxies. We can com- 
pare this result with similar studies of SO galaxies base d 
on: (1) optical and NIR photometry l|Peletier et al.lll999l ). 
and (2) integra l field unit (IFU) optical spectroscopy 
l|Sil'chenkdl2006h , as well as comparing the ages and metal- 
licities of specific galaxies in our sample with spectroscop- 
ically derived values. Since the central ages of early-type 
galaxies, which may in clude SOs, are typica l ly older for 
more massive galaxies llCaldwell et ail 120031; I Nelan et alj 



20051; ISanchez-Blazquez et al l 120061; iGraves et all l20Qg[ : 



Isanchez-Blazguez et all |2007| ;~ Smith et all 120091 ) ~we must 
compare galaxies of similar stellar mass (absolute luminos- 
ity). The typical stellar masses of our sample galaxies are 
similar to those of the Peletier and Sil'chenko samples. Thus, 
stellar masses cannot be responsible for any age differences 
between the samples. 

From broadba nd optical and NIR photometry, 
iPeletier et all (|l999l ) found a mean bulge age of 9 ± 2 Gyr 
(at 1 bulge effective radius) for a sample of twenty SO and 
early spiral galaxies (to Sbc). The "bulges" of lPeletier et al.l 
(1999) are closest in observed radius to our 2nd binning 
region, for which we find a mea n age of 4.6 ± 2.4 Gyr. 
While the IPeletier et al. I l|l999h ages are substantially 
greater than ours, these authors did not include the seven 
bluest, youngest bulges in their sample, three of which are 
SO galaxies. As well, they studied highly inclined galaxies 
and masked out the disk from their extracted ages and 
metallicities, hence their results truly apply to bulge light 
at 1 R e , while ours apply to bulge and disk light mixed 
together within 0.5 to 1.5 R e . Finally, the difference in 
mean ages may be partly attrib uted to t he difference in 
SP m odels. IPeletier et ail (|l999h used the Vaz dekis et alj 
i|l996l ) SSP models but note that the IWorthevi (|l994T ) 
models give a significantly smaller mean age (2 Gyr) for 
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their sample. The latter would be in closer agreement with 
our result. 

ISil'chenkol (|2006l ) compiled an IFU spectroscopic sample 
of 58 lenticular galaxies and used the Thomas et al (2003) 
models for age estimates. She reports median ages within 
the unresolved nuclei of 3.7 and 6 Gyr for galaxies in sparse 
and dense environments, respectively, while the median ages 
for 'bulges' (extracted from an annulus of 4 to 7 arcsec) are 
4.8 and 8.3 Gyr for sparse and dense environments. With the 
exception of the few Virgo galaxies, our sample shares the 
sparse environment of Sil'chenko's galaxies, and our central 
ages are extracted from a region between the unresolved 
nucleus and the 4-7 arcsec b ulge. Our a ge di stribution is 
thus a closer match to that of ISirchenkol 1 20061 ). 

We also wish to perform a comparison of metallicities, 
for which the mean and median central values of our sample 
galaxies are [Z/H]= 0.2 ± 0.3 dex and 0.3 dex, respectively. 
The median value is preferable here, as it is more robust 
to the limit we place on me tallicities a t +0.5 dex. Our re- 
sults are in agreement with ISirchenkol (|2006l ), who finds a 
median metallicity for her nuclei and bulges of +0.4 dex 
and +0.2 dex , respe ctively. Visual inspection of Fig. 3a in 
iPeletier" et alJ (|l999h reveals a mean metallicity near solar, 
consistent with our mean value of [Z/H]=+0.0 ± 0.3 dex for 
our second radial bin. 

Although our sa mple does not o verlap with either 
IPeletier et all (|l999l ) or lSil'chenkol (|2006h ■ independent ages 
and metallicities from spectroscopic studi es offer direct com- 
pariso ns for a few galaxies in our sample. iKuntschner et all 
(2006) extracted Lick index line strength measurements for 
UGC 5503. Co mparison o f their H/3 and Fe5015 index mea- 
surements to a IVazdekij (| 19991 ) stellar population model 
yields an age of 1.6 Gyr for the inner 25". Our average ex- 
tracted age for this region (compara ble to our first two ra - 
dial bins) is 1.3 Gyr. For UGC 10048. ICaldwell et all (|2003h 
estimated an age range of 4.9 to 7.9 Gyr, depending on 
the indices used. They used a 3" wide long slit spectrum 
with variance weighting, thus it is difficult to compare spe- 
cific regions in the galaxy. Nevertheless, our extracted age 
is slightly younger, ranging from 4.5 Gyr in the inner re- 
gion to 3.7 Gy r in th e outer region. The metallicity from 
ICaldwell et al] l|2003h ranges from -0.1 to 0.1 while ours 
ranges from 0.1 in the inner region to -0.3 in the outer region. 
ICaldwell et ail (|2003h also studied VCC 1614, for which our 
extracted age for this galaxy ranges from 1.4 Gyr in the in- 
ner region to 5.5 Gyr in the outer region. This is consistent 
with their age of ~1.8 Gyr. Our metallicity ranges from -0.4 
in the inner region to -1.2 in th e outer region, consis tent 
with the metallicity of ~ -0.4 for ICaldwell et all l|2003h . In 
general, our extracted ages and metallicities are similar to 
those derived from spectroscopic studies. 

Fig. [12] shows the center of several galaxies lying 
off the model grid toward larger r-H values. The col- 
ors of these off-grid points may result from dust redden- 
ing. In man y SO galaxies, dust affects mainly the cen- 
tral regions JPeletier et all Il999l ; IWikl nd fc Henkell l200ll ; 
iFritze v. Alvenslebenl 120*04 ). although dust lanes are ob- 
served in 3-color optical images in the outer regions of a few 
galaxies in our sample. While our galaxy colors are corrected 
for Galactic foreground extinction, internal extinction is ig- 
nored. To estimate the effect of internal extinction, a redden- 
ing vector for a foreground screen dust model with A„ = 0.3 



is shown in the upper left corner of the color-color diagrams. 
Thus, if unaccounted for, intrinsic extinction biases our ex- 
tracted ages and metallicities towards highe r values. This 
effect is clearly demonstrated and discussed in lPeletier et al.l 
(1999). Their Figure 3a shows a large shift in colors, parallel 
to the reddening vectors, from the center of the galaxy to one 
R e ; ages could shift from ~ 11 Gyr to ~ 9 Gyr and m etallic- 
ity from [Fe/H] > 0.5 to ~ 0.0. IPeletier et al] 0999) found 
a dust signature in HST images in almost all SO galaxies in 
their sample. 



7.2 Radial Age Trends 

While the inner regions of our SO galaxies have light- 
weighted ages of 3.7 ±1.8 Gyr (the error is rms scatter), 
the outer regions are remarkably heterogeneous. For those 
we find a mean age of 6.1 ± 4.0 Gyr. The ± 4.0 Gyr scat- 
ter in age exceeds the observational error. The larger age 
spread in the outer regions is notable by comparing the top 
and bottom panels of Fig. 1141 which shows all galaxies in 
our sample (as mentioned before, we use the H-band data 
for the 11 overlapping galaxies). Many galaxies have outer 
ages that exceed our limit of 13.8 Gyr. The increased hetero- 
geneity in outer region ages compared to the inner regions 
naturally creates heterogeneous radial age trends, as illus- 
trated in Fig. [T2j We find an increase in age with radius 
from the inner to outer radial regions for 58% of our sample 
of 59 galaxies, a decrease for 19%, and a change of less than 
1 Gyr for 23%. In Fig. 1141 there is a clear tail in the outer 
region for high ages. As expected, this tail is also apparent 
in the distribution of age differences from the inner to the 
outer region, shown in the top panel of Fig. [15] Excluding 
this tail, many galaxies have similar ages between the inner 
and outer regions. 

We find that the outer regions for 24% of our galaxies 
show a substantial increase in light-weighted age from the 
center of the galaxy outward, resulting in outer regions older 
than 10 Gyr. We refer to these galaxies as "old outskirts" 
galaxies. Specifically, we define an "old outskirts" galaxy to 
be one for which the integrated age of the "outer" region 
(i.e., all radii beyond 1.2 R e ) is greater than 10 Gyr. Note 
that some galaxies with ages < 10 Gyr in one or more outer 
radial bins are not classified as "old outskirts" galaxies be- 
cause the integrated age for all radii beyond 1.2 R e does not 
exceed 10 Gyr. However, we include two galaxies as "old 
outkirts" whose H-band derived ages in the outer region 
are just under 10 Gyr, but the J-band derived outer region 
ages exceeds 10 Gyr (see t|7.2.ip . Examination of optical 
color images for those galaxies shows no prominent features 
(i.e. bars, nascent spiral structure or rings) in their outer 
regions, except for UGC 9713 (bar) and UGC 10112 (faint 
ring). Color images for galaxies with old outskirts are shown 
in the electronic version and are indicated by an "O" above 
the image. 

We have veried that the aforementioned results are ro- 
bust to the definition of the outer region by varying it from 
1.2 to 3.5 R e . We see no change and, thus, conclude that the 
definition of the outer cutoff does not change strongly affect 
our results. 

We can compare our age trends with those from simi- 
lar studies. Early photometric studies of SO galaxies found 
mostly negative age gradients, but covered only the in- 
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ner disk and used different color com b inations than ours 
Caldwelj[l983l : iBothun fc GreggHl990T l. iPeletier fc Balcellsl 



19961 ) find an age difference of less than 3 Gyr between the 
bulge and inner disk for most galaxies in their sample, sim- 
ilar to our peak at zero age difference. While they did not 
have many old outskirts galaxies, their study included early- 
type galaxies up to Sb, with only 8 pure SO galaxies, two of 
which have older ages in the outer regions. Detailed studies 
of individual SO ga laxies have found both negative and pos i- 
tive age gradients (TTikhonov et alj|2003l ; iNorris et aLll2006h . 
Finally. ISirchenkoT i 2006T ) found that, independent of envi- 
ronment, SO galaxies contain a nuclear region that is younger 
than the bulge region. We show here that this trend contin- 
ues outward for some galaxies in our sample, resulting in a 
large age gradient for those. 

We can comp are our age gradients with those of 



U (I1996J) . also discussed in t )7.1l and those of 

u die 



Fisher et al. 

Rawle et al.l (|2010T l. who studied a sample of 19 early-type 
galaxies, including 11 SOs, in the Shapley Supercluster core. 
For our galaxies, we find a mean Aagc/Alog(r/R e ) = 2.3 
± 4.6 Gyr dex -1 , with a media n of 0.6 Gyr dex -1 . For a 
sample of 9 edge-on SO galaxies, iFisher et al.l (| 19961 ) found 
Aage/Alog(r/R e ) = 5.0 to 7.0 Gyr dex -1 (depending on 
assumptions) for the bulge and Aage/Alog(r/R e ) < 1 Gyr 
dex -1 for the inner disk. Even though we sample light from 
both the bulge and disk components due to the face-on na- 
ture of ou r sample, we do fin d that our age gradients are 
consistent. iRawle et al] (|2010h found a lower mean age gra- 
dient of Aage/Alog(r/R e ) =-0.10 ± 0.04 Gyr dex" 1 , but also 
point out the existence of several individual galaxies with 
significant gradients in either direction, as we have found. 



7.2.1 Caveats for Old Outskirts SO Galaxies 

A striking finding of this study is the significantly old ages in 
the outer regions of 14 SO galaxies. We discuss these galaxies 
in more detail below along with potential caveats, whether 
observational or model dependent, for our inference of old 
outer regions. To assist our discussion, we show in Fig. \W\ 
color-color diagrams for these 14 galaxies. 

We now test for observational consistency between J 
and H band data for the old outskirts galaxies. Based on 
H-band data, three galaxies, UGC 10163, UGC 10391 and 
UGC 9280, would not qualify as "old outskirts", but do so 
by virtue of J-band data. The H-band color-color diagrams 
for UGC 10391 and UGC 9280 show that those galaxies do 
have old outskirts, such as ages > 10 Gyr in at least 2 outer 
radial bins (but narrowly fail to be classified as "old out- 
skirts" galaxies based on their integrated ages beyond 1.2 
R e ). Examination of the color image of UGC 10163 (see the 
electronic version) reveals dust in the outer region of the 
galaxy which artificially bias towards old ages. To be con- 
servative, we eliminate this galaxy from our sample of old 
outskirts galaxies. Two galaxies in our sample with J and H 
band data, UGC 9713 and UGC 10112, qualify as an old out- 
skirts galaxy using ages extracted from both bands. Thus, 
our final sample of old outskirts galaxies includes four SOs 
with data in both the J and H bands: UGC 9280, UGC 9713, 
UGC 10112, and UGC 10391. The color-color diagrams for 
those galaxies (Fig. I16[) are very similar between J and H 
(the difference in UGC 9713 in the last bin is due to a stricter 
cutoff in the J-band). 



Not only do old outskirts galaxies have old outer re- 
gions, they also have a positive age gradient. Because of the 
tilt of the model grid in Fig. 1161 the outer region ages and 
age gradients are highly dependent on the r-H color, even 
though the g-r color is more sensitive to age than r-H. Thus, 
our results hinge on the observed large decrease in r-H color. 
This underscores the importance of NIR colors in extracting 
ages from photometry. 

A major concern for the identification of faint outer re- 
gions in galaxies is the sky subtraction, which, as discussed 
in £14.41 is more problematic in the NIR than the optical. 
Possible systematic sky errors include a variation of sky 
level across the image field, a narrow field of view for the 
galaxy size (never a problem here), and a bias in the sky 
level assessment. We discussed in H4.4I NIR systematic error 
estimates based on the variation in sky level from four boxes 
along the perimeter of the image, which is included in the 
error bars in Fig. 1161 Provided our sky error estimates are 
correct, altering the sky background by a systematic error 
in most of these galaxies would not yield an age below 10 
Gyr; that is, old outer regions remain. Finally, any error in 
our methodology in sky background estimates should affect 
all of our galaxies equally. 

It is also unlikely that the presence of dust has skewed 
our results of old outer regions. The reddening line shown in 
the top left of Fig. ll6l shows the direction for dust reddening, 
which would make both g-r and r-H redder. Thus, a large 
amount of dust reddening would be necessary to explain 
the observed ol d ages. While SO ga laxies likely suffer some 
dust extinction. IPeletier et al l (|1999| ) find from HST images 
that the light should be free from dust for most SO galaxies 
beyond r > 1 R e . Examination of SDSS color images of 
old outskirts galaxies shows no visible dust lanes or nascent 
spiral structure in the outer regions of these galaxies. 

We have also examined photometric zeropoint errors. 
However, because this would affect all radii equally, no 
amount of shift in any direction would negate the trend of 
increasing ages with radius for old outskirts galaxies. 

Our derived ages are notoriously subject to modeling 
uncertainties, especially at NIR fluxes that are dominated 
by il l- constrained late evolutionary phases dEminian et al.l 
120081 ; ICharlot fc Bruzuall l20ld ; IMacArthur et al.l |2010| 1.A 
major challenge for SP models is the treatment of the 
TP-AGB phase which affects mostly low ages and high 
metallicities, but old ages are unchanged. However, col- 
ors for the outer regions of old outskirts galaxies fall off 
of all model grids, suggesting an inconsistency in obser- 
vations and the models. While some of the offset may be 
due to sky subtraction error, this cannot explain the off- 
set entirely. Thus it is possible that zero point errors in 
the SP models are at play and may affect our results. 
Models have been known to predict excessive galaxy ages 
dVazdekis et al.ll200lj; IMacArthur et al.|[20o3 : iThomas et all 
120051 ; iMoorthv fc HoltzmanlbOOd l. Age dating of old stellar 
populations is influenced by uncertainties in the T e // and 
[Fe/H] scale of giant branch stars (a significant fraction of 
the continuum light of old stellar populations, even in the 
blue, is provided by giant stars) as w ell as the luminosity 
function of the upper red giant branch (|Schiavon et al ] |2002l . 
hereafter S02). For example, S02 found that the unrealisti- 
cally large derived age for the globular cluster 47 Tucanae 
from comparison of its observed integrated spectrum to pop- 
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ulation synthesis models may be partially explained by the 
discrepancy between the observed luminosity function of the 
upper red giant branch of the cluster and the lower expected 
number based on synthesis models. Indeed, the model pre- 
dictions in S02 for (B-V) are too blue compared to the ob- 
servations of 47 Tucanae, so that an older age is needed to 
reconcile with the observed color. Despite these concerns for 
the model-dependent derivation of old ages, it is unlikely 
that the result of old outer regions in a sub-sample of SO 
galaxies will be revised with improved models. 

Finally, our reliance on SSP models to extract ages 
should not affect this result. Since young stars dominate 
the light, an SSP age greater than 10 Gyr indicates the lack 
of any significant star formation for quite some time. There- 
fore, no star formation history is necessary. 

Given all of the above considerations, we feel that the 
result of old galaxy outskirts is robust in most cases. 

7.3 Age Trends With Galaxy Properties 

We have established that SO galaxies exhibit a variety of 
radial trends in age, with a substantial population having 
large positive age gradients. We now inquire whether these 
age trends, as well as the mean ages, are correlated with 
global properties such as total mass and light concentra- 
tion. To that end, we have separated our sample of 59 SO 
galaxies into two roughly equal groups of low and high mass 
galaxies, straddling lxlO 11 Mq. The colors in each radial 
bin are averaged for all galaxies in each group and are plot- 
ted in the g-r vs r-H diagram (left panel) and g-r vs r-J 
diagram (right panel) in Fig. 1171 The solid colored error bars 
are the standard error in the mean for each radial bin. The 
black dotted error bars in the first radial bin denote the 
la scatter in color for galaxies in that bin. In both panels a 
separation in ages is apparent between the two mass groups, 
such that the stellar populations of less massive galaxies are 
younger on average than those of more massive galaxies at 
all radial bins; this effect is more pronounced in the J-band 
sample. Although the sample is not statistically complete, 
ac counting for pos sible environment-age trends (as observed 
m ISiPchenkdlioblih would only strengthen our mass-age re- 
lation. A KS two-sample test was applied to the high and 
low mass age distributions for all radial bins combined for 
all galaxies in our sample. The hypothesis that the two sam- 
ples are drawn from the same population in mean age can 
be rejected, with a p-value of 1.7x10 , confirming what is 
apparent in Fig. 1171 The results of this and other KS tests 
are summarized in Table 7. When only the inner regions are 
considered, the ages of the two mass groups can still be re- 
jected as a single population ( at the 97% signif icance level). 
This result is consistent with ISil'chenkd (120061 ). who found 
that the stellar populations of more massive SO bulges are 
older than those of less massive bulges. However, the bias in 
mass and environment of our sample (see [JSJ may weaken 
age-mass relations. 

We have also sorted galaxies into high and low concen- 
tration groups, using a separation at the median value of 
the concentration index, C28 = 4.7. As above, the ages at 
each radial bin were combined for each galaxy in the two 
concentration subsamples. Fig. [TH] shows a color-color dia- 
gram for the J- and H-band samples with the average colors 
of each concentration group averaged. A separation in age is 



evident, such that the stellar populations of centrally con- 
centrated galaxies (C28 > 4.7) are older than those of less 
centrally concentrated galaxies. As with mass, the J-band 
sample, which does not include the Virgo galaxies, shows a 
clearer separation in age. The hypothesis of a similar parent 
distribution in age between the two concentration groups 
can be rejected, with a p-value of 1.9xl0 -4 , when all radial 
bins are combined. We do not find a statistically signifi- 
cant (within 1 a) age difference in the outer regions alone 
for different concentration groups. To show the mass-age 
and concentration-age relations without the constriction of 
groups, we plot in Fig. 1 191 age versus mass and concentration 
for the inner and outer radial regions. 

Because low concentration galaxies tend to be less mas- 
sive than higher concentration ones, as seen in Fig. [5J the 
age trends found for both mass and concentration could be 
intertwined, with only one of the variables being the driving 
factor. To test which of mass or concentration is more funda- 
mental, we further subdivide our high mass sample into low 
and high concentration subsamples, and do the same for the 
low mass sample. In Fig.J^njcolor-color diagrams are plotted 
with galaxies separated by mass as in Fig. 1171 however this 
time only high concentration galaxies are plotted in the left 
panel and only low concentration galaxies are plotted in the 
right panel. Fig. 1211 reverses the analysis. Figs. 1201 - 1211 show 
a larger separation in mean age for the low concentration 
and low mass galaxies. While this may suggest that both 
parameters drive the correlation in age, these eyeball tests 
are not necessarily statistically sound (at least within 1 a). 
Combining mass and concentration, we find that galaxies 
with both low mass and low concentration are younger in 
mean age than the rest of our sample in the inner region 
(2.9 Gyr versus 4.2 Gyr). 

Our sample of SO galaxies is a morphological mix, 
with some galaxies showing more spiral features than oth- 
ers. We now examine how radial age trends are depen- 
dent on morphology. The distributions of age differences for 
non-transition galaxies (featureless disk) and for transition 
galaxies are plotted in Fig. [TS] in the middle and bottom 
panels respectively. There is a positive age gradient for fea- 
tureless disk SOs which is clearly stronger than that for tran- 
sition galaxies. Transition galaxies do not appear to have the 
same tail of large age differences as other galaxies. This is ex- 
pected if transition galaxies have a growing spiral structure 
where young stellar populations are located. A KS test for 
the age differences of transition and non-transition groups 
reveals statistically different populations at a greater than 
3a confidence level. 

We note that old outksirts galaxies do not exhibit a 
preference for any particular mass or environment, but may 
for concentration. Also, age gradients, do not demonstrate a 
correlation with mass or concentration. 

7.4 Metallicity Trends 

Negative metalli city gradients have been previously found 
for SO galax i es dFisher et all [l996l: iTamura fc Ohtal 120031 : 
iRickes et alJ l2009i: iRawle et al.l 12010) . From their spec- 
troscopic sample iFisher et al.l ([1996) found a mean of 
A[Fe/H]/Alog(r/R e ) = -0.9 to -0.7 (depending on assump- 
tions) for their bulges a nd A[Fe/H1/A(r fcpc ) = -0.04 to - 
0.06 kpc~ for the disks. IRawle et all (|2010) found a mean 
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metallicity gradient of A[Z/H]/Alog(r/R e ) = -0.13 ± 0.04. 
We carry our measurements typically 2.5 to 5 times further 
in radius and look for correlations with global properties. 
Since colors alone do not yield information about non-solar 
element abundance ratios, we examine the trends of mean 
metallicity in SO galaxies based on the scaled solar modeling 
described in §6. 

Fig- HHshows that virtually all SOs have significant neg- 
ative metallicity gradients. We find a mean metallicity gra- 
dient A[Z/H]/Alog(r) = -0.6 ± 0.5. In terms of physical 
radius, the gradients have a mean A[Z / H]/ A(rk pc ) = -0.1 
± 0.3, consistent within the errors with the above studies. 
Recall that the central-most binning region is avoided in the 
calculation of gradients. This largely alleviates any artificial 
metallicity bias due to dust, if any, in galaxy centers. 

The correlation of metallicity with mass and concen- 
tration can be assessed in the same way as for age, i.e., by 
dividing the sample into high and low mass and concen- 
tration groups. For the trends of inner region metallicities, 
eleven galaxies lie off the model grid and are thus set to the 
maximum metallicity of 0.5 dex; this would most likely have 
a weakening effect on metallicity trends. As well, accounting 
for the fact that all of the observed lowest mass and lowest 
concentration galaxies are in the Virgo cluster (high density 
environment) would only strengthen the mass and concen- 
tration trends that we observe, because galaxies in higher 
density environme nts are found to ha ve higher metallicity 
(at a given mass) (jCooper et al.ll2008l ). 

Examination of Figs. [17] and [18] shows a strong sepa- 
ration in metallicity for the high versus low mass and con- 
centration groups such that low mass or low concentration 
galaxies have lower metallicities. High and low mass galaxies 
have statistically different metallicity distributions in both 
the inner and outer regions (KS test reveals a p-value of 
6.77xl0~ 5 for the inner regions and 0.01 for the outer re- 
gions). Studies of mass-metallicit y relations in SO galax- 
ies are abundant dGarnettl |2002|; iKauffmann et al.l 20031; 



MacArthur et alj|2004 iTremonti et al.ll2004l ; iThomas et all 
20051 ). but most of these studies report a central metallic- 
ity (< 1 disk scale length) or a global metallicity (such as 
studies at high z where the aperture encompasses the entire 
galaxy) so that a comprehensive mass-metallicity relation in 
the outer parts of SO galaxies is still missing. The metallicity- 
concentration trends we see may be related to the increase 
in metallicity with early Hubble type (M04). Metallicity 
gradients, on the other hand, do not demonstrate a correla- 
tion with mass or concentration. The mass-metallicity and 
concentration-metallicity relations are also demonstrated in 
Fig. [22] for the inner and outer radial regions. 

We try again to separate the effects of mass and con- 
centration on metallicity. Examination of Figs. [5IJ] and 1211 
shows a significantly larger separation in mean metallicity 
for the low concentration and low mass galaxies. For the 
inner regions, the hypothesis that high and low mass galax- 
ies of low concentration are from the same population and 
that high and low concentration galaxies of low mass are 
from the same population can be rejected, with a p-value 
of 2.5 x 10~ 4 for both. There is no statistical difference for 
the high mass galaxies. That the high mass sample does 
not show a difference in metallicity for concentration groups 
may either mean that concentration does not play a large 
role in metallicity trends or that it does not play a role in 



the range of concentration covered (C28 > 3.8). Regardless 
of the driving force for the metallicity relations, we find that 
galaxies with both low mass and low concentration have a 
lower mean light-weighted metallicity in the inner regions 
than other galaxies in our sample. 



8 DISCUSSION 

We have examined the connection between the globally av- 
eraged light-weighted ages and metallicities of SO galax- 
ies and other global properties, such as stellar mass and 
concentration. We have also investigated age and metallic- 
ity gradients in SOs. We find that at all radii, SO galax- 
ies with lower mass and lower concentration have younger 
light- weighted ages and lower light- weighted metallicities. 
We have also found that virtually all SO galaxies in our 
sample have negative metallicity gradients, with an average 
value of A[Z/H]/Alog(r/R e ) = -0.6 ± 0.5, consistent with 
published results. On the other hand, the radial behavior 
of age in SOs is heterogeneous; both positive and negative 
age gradients are found. We observe an increase in age with 
radius for 58% of our sample, a decrease for 19%, and lit- 
tle change (< 1 Gyr) for 23%, and a mean age gradient 
of Aage/Alog(r/R e ) = 2.3 ± 4.6 Gyr dex' 1 . For 24% of 
our sample, we find populations with substantially old light- 
weighted ages (> 10 Gyr) in the outer region of the galaxy. 
We now consider the implications of these results. 



8.1 Mass and Concentration Relations 

The mass-metallicity trends that we observe extend into the 
outer regions of the galaxy, which may be dominated by 
a disk component, but we can not be certain without full 
bulge-disk decompositions of our galaxies. We assume here 
that the inner regions are bulge dominated, but this, as well, 
is not certain without decomposition analysis. 

Both high galaxy concentration and high galaxy mass 
point to a higher bulge mass. We speculate whether a rela- 
tion between the bulge mass with age and metallicity in the 
inner and outer regions favor a particular model for SO for- 
mation. Observations consistent with bulge growth through 
internal (i.e. driven by disk instabilities) and external (i.e. 
driven by satellites) secular evolution hav e been noted in 
vario u s studies of lat e type spiral galax ies (Courtcau et al 



1996 1; ICombel |2000|; lEllis et al.l pOOll; iKannappan et al 



2004 ; iKormendv fc KennicuttJ 2004l ~ Barwav et al. I l2007t ). A 
small bulge, such as in late type spirals or low mass, low 
concentration SOs, is more likely to contain associated star 
formation from secul ar evolution, eve n if most of the mass 
has an earlier origin (Ellis et al. 2 00lf ). In deed, the spectro- 
scopic analysis of stellar populations by iMacArthur et ail 
(2009) reveals that ~ 20% of the mass of late type spiral 
bulges consists of stars less than lGyr (although, 70% of 
the bulge light comes from these young stars) . 

The merging scenario provides another explanation for 
the age and metallicity trend with mass and concentration in 
the inner region. It has been shown analytically that positive 
correlations of mass with average age can be a natural con- 
sequence of hierarchical merging |Neistein et alj|2006h . As 
well, detailed semi-analytic modeling of hierarchical merg- 
ing that includes enhanced feedback processes such as AGN, 
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yield more extended star formation histories in less massive 
galaxies, leading to posi tive correlations of mass with aver- 
age age and metall i citv (|Cole et al.ll 19941 : iBower et al.ll2006l : 
ICroton et al]|2006l ; iDe Lucia et al.ll2006l ) . 

Thus, the fact that we find mass-metallicity relations 
for the inner regions of SOs does not discriminate between 
models where the bulge is formed through major mergers 
versus secular evolution models; they both predict age-mass 
and mass-metallicity relations for the bulges of SOs. Note 
that such relations also naturally result in a "monolithic 
collapse" scenario. 

However, we see mass-metallicity relations in the outer 
regions as well. Can this point toward a part icular forma- 
tion mechanism? I Yoachim fc DalcantorJ (|2006h propose that 
increased outflows in low mass disks at early times (in sub- 
units before assembling into a disk) allows a mass-metallicity 
relation in disks to be established early so that hierarchi- 
cal merging may be able to explain the observed mass- 
metallicity relation in the outer regions of SO galaxies. It 
is unclear whether such a relation can be explained by sec- 
ular evolution considering star formation must be efficient 
enough in the outer regions for enrichment to occur. While 
the mass-metallicity relation observed in the outer regions of 
SO galaxies may be a useful discriminant of formation mech- 
anisms, predictions from both merger and secular evolution 
models have yet to consider this scenario. 

The parallel between our finding of younger ages for low 
mass and low concentration SOs and a similar finding that 
late type spirals are younger than earlier types (M04) is con- 
sistent with a similar formation picture for SOs and spirals. 
We have demonstrated in i]7.3l that the structural transition 
in "transition" galaxies correlates with a transition in SP 
properties between SOs and spirals; as with the spiral galax- 
ies observed by M04, transition SO galaxies do not show 
the tail of increasing ages that we observe for non-transition 
galaxies. In fact, the majority of spiral galaxies are shown to 
have decreasing ages with radius (M04). The transition SO 
galaxies thus appear to bridge the gap between SO and spi- 
ral galaxies, supporting a formation connection between the 
two groups. Transition galaxies are likely candidates for de- 
scendants of spiral galaxies that have lost or exhausted their 
gas. However, the fact that many transition galaxies in our 
sample have large stellar masses is incon sistent with the pro- 
posed mass dichotomy in SO formation (Jorgcnse n fc Franxl 
1 19941 : iMehlert et"aDl2003l : iBarwav et"ai1l2007T ). namely that 
only low mass SOs originated from spirals, while the high 
mass SOs have a formation history similar to ellipticals. 

8.2 Old Outskirts Galaxies 

We have found a subsample of galaxies with large positive 
age gradients and significantly old outer regions (> 10 Gyr). 
These galaxies appear to form a separate evolutionary class 
from the rest of the sample, showing no preference for mass, 
concentration, or environment. Their physical appearance 
is largely featureless. We have discussed potential observa- 
tional and model dependent caveats for the result of old out- 
skirts in i|7.2.1l Assuming this result is robust, we explore 
tentative implications of galaxy evolution. 

Until we gather more information from light decom- 
positions, we can only speculate on the nature of these 
old outskirts galaxies. The old populations in the outer re- 



gions either belong to a disk or a spheroid component. In 
the case of a spheroid membership, the galaxy may have 
a dominant spheroid population with an embedded disk 
so that the radial increase in ages is due to the change 
in luminosity-weighted ag es from a disk-dom inated to a 
spheroid-dominated region. iPeletier et ail (|2007l ) found that 
the minor axis, and thus the spheroid component, of 24 
early-type (Sa) galaxies has uniformly old SPs. It is this 
old spheroidal component that might be dominating in the 
outer regions of old outskirts SOs. In this scenario, old out- 
skirts galaxies would be closely tied to elliptical galaxies. In 
fact, the latter are kno wn to harbor small embedded disks 
|Krainovic et all 120081 . M10), hence the old outskirts SOs 
would tend to straddle between ellipticals and SO galax- 
ies. Luminous r ed halos have also b e en observed around 
spiral galaxies jLequeux et all Il998l: IZibetti et al l 120041 : 
IZackrisson et al.l 120061 ). However, IZackrisson et alj (2006) 
found that old populations cannot explain the red colors of 
many spiral halos without invoking an extremely bottom- 
heavy IMF. 

If, instead, the outer regions of old outskirts galaxies 
are dominated by a disk component, their stars would have 
formed long ago. Given uncertainties in SP modeling, set- 
ting an absolute age on the formation of the disk is difficult, 
but in old outskirts galaxies the stars dominating the light 
in the outer regions clearly formed long ago (possibly > 10 
Gyr). Moderately old disks (~ 7 Gyr) have been predicted 
in theoretical models where t he main infall phase prec edes 
the onset of star form ation (|Ferguson fc Clarke! 1200 li ). As 
well, iKoo et all (|2005l ) find several very red luminous disks 
at high redshift, implying that at z = 1 not all massive disks 
are young; some old, massive SOs may have already existed 
in the field. For reference, stars formed at z = 1 would be 
~ 7 Gyr old at the present epoch. The presence of a disk 
component implies that no major, disk-destroying event has 
occurred since the formation of the dynamically fragile stel- 
lar disk; stellar disks are thought t o be destroyed in merger of 
galaxies with a mass ratio > 1/4 jMihos fc Hernquistlll996l : 
ISteinmetz fc Navarrol [20021 ). Thus, the very old ages of the 
galaxies' outskirts (possibly > 10 Gyr) place a constraint on 
models of hierarchical merging if the outskirts have a disky 
origin. 

Previous studies of spirals with low present day B/D 
suggest that hierarchical models may req uire no major 
merg er since the formation of the stellar disk (|Weinzirl et al.l 
2009). In principle, low B/D spirals may be produced 
through gas rich mergers, which produce a new disk in th e 
merger remnant (|Hopk ins et al. I l2009l : iRichard et aT]|2010h . 
However, our observations of old outskirts galaxies cannot 
be explained through these mergers because new stars will 
necessarily form in the disk and would then dominate the 
light. Thus, the requirement of no major merger since z ~ 2 
(using an age of 10 Gyr), if the outer regions are indeed a 
disk, is necessary to explain our observations. 

Alternately, old stars may have somehow been redis- 
tributed into the outer regions. However, because the inte- 
grated light is heavily weighted toward young SPs, the old 
light-weighted ages in the galaxies' outskirts imply that vir- 
tually no young stars exist in this region. An accretion of 
old satellites to the oute r disk may be possible. For exam- 
ple, Ide Jong et al.l (|2007T l found that a stream around M83 
is made of an old stellar population. Additionally, a migra- 
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tion of ol d stars from the cen ter outward is predicted in the 
models of iRoskar et al.l l|2008l ) through spiral corotation res- 
onance scattering, although the lack of spiral structure in SO 
galaxies may make this difficult. 

The large ages and low metallicities in the outer re- 
gions of old outskirts galaxies rese mble those of thick disks 
observed around edge-on galaxies l|Dalcanton fc Bernstein! 
120021; lYoachim fc Dalcantonl I2006T ) and in the Milky Way 
i|Gilmore et al.lll989l ). Perhaps the populations in the old 
outskirts of galaxies belong to a thick disk. Indeed thick 
disks may form during a gas-rich merger event and an em- 
bedded thin disk forms afterward (jRichard et al.ll2010r ). If a 
galaxy has its gas stripped shortly after this merger event, 
perhaps part or all of the young thin disk will be kept from 
forming and the galaxy will contain old outskirts. Unfor- 
tunately, our galaxies are seen face-on so it is difficult to 
determine whether the outer populat i ons be long to a thick 
disk. We note that the lRichard et all l|2010l 'l models invari- 
ably show a major starburst occurring during the merger. 
Thus, mergers like these must be old if they are to explain 
the outer populations of old outskirts galaxies. 

Most disk/spiral galaxies s how a decrea se in 

light- weighted a g e wi t h radius feell fc de Jonel |2000| ; 
iMacArthur et all 12004 iRoediger et alj l201Ch . " a trend 
opposite to what is observed in old outskirts galaxies. 
The positive age gradients of old outskirts galaxies, which 
may be due to a difference in bulge and disk age and/or 
to an age gradient in the disk, can offer clues to their 
formation. If the bulge (or the population dominating the 
light in the bulge) formed after the disk, some mass in the 
bulge would need to be created without either destroying 
the fragile disk or adding star formation to it. Internal 
secular evolution in gaseous galaxies is expected to fuel star 
formation in a centralized disk. Internal secular evolution 
may play a role in creating the young bulges in old outskirts 
galaxies, but other processes would be required to explain 
the suppression of s tar fo rmation in the disk. Alternately, 
iHernauist fc Mihosl j 19951 ) show that even a 10-to-l merger 
can drive up to 50% of the primary galaxy's gas into its 
center. Not only will this help to build up a young bulge, 
but it will suppress star formation in the disk. 

A positive age gradient in the disk component itself 
could be due to either the original formation of the disk 
or to a later event that transformed the galaxy into an 
SO. Proposed theor ies of SO formation through external gas 
removal processes (IQuilis et al. I l200d ; iKawata fc Mulchaevl 
120081 ; iKronberger et a!T 20081 ) call for an outside-in strip- 
ping of the galaxy. The resultant galaxy will have a lower 
light-weighted age inside the truncation radius than beyond 
it. Additionally, enhancement of star formation from sev- 
eral gas removal processes is expected in the central regions, 
again yielding an increase in age with radius. The events 
that are expected to bring about this age difference, how- 
ever, are also ex pected to be short-live d with a timescale 
around 500 Myr IjKronberger et al.ll2008l ). To create the ob- 
served large positive age gradient, star formation must be 
sustained in the inner regions for a very long time, possibly 
through minor accretion. The age increase in old outskirts 
galaxies' disks may instead be due to the original formation 
of the disk. 

Posi tive age grad i ents in disks are not exclusive to SO 
galaxies. I Taylor et al.l (|2005l ) found that the outer regions 



of some late type spiral galaxies are redder and the inner 
regions are bluer than other galaxies in his sample, suggest- 
ing that this could be an indicator of outside-in formation, 
with significant recent star formation in the inner regions 
and minimal recent star formation in the outer regions. A 
positive age gradient is predicted in simulations of the dissi- 
pational collapse of gas embedded in a s pherical dark matte r 
halo through scattering of disk stars ([Roskar et al.| [2008). 
Beyond a nominal radius, R oskar et al.l (|2008T l find an in- 
crease in mean stellar age for their model disk galaxy. Spiral 
arms at some point in the history of the host galaxy are a 
necessary component for this model. 

The scenario of the thick disk forming first may also be 
used to explain radial trends in old outskirts galaxies. Most 
likely, several processes are involved, i.e. perhaps the forma- 
tion of the thick disk is halted through a partial stripping 
mechanism. Young disks embedded in the thi ck old disk and 
the ou tside-in disk growth in the models of iRichard et al.l 
(2010) may explain the trend in radius toward older ages 
observed in old outskirts galaxies. 

In summary, the observed mass-metallicity trends in SO 
galaxies are consistent with a merging scenario for the for- 
mation of the bulge, but the cause of the relation that we 
observe in the outer regions has yet to be considered in mod- 
els. Our results, in particular those of transition galaxies that 
straddle SO and spiral galaxies, are consistent with the for- 
mation of SO galaxies through the transformation of spiral 
galaxy disks. It is possible that old outskirts galaxies require 
no major merger to have occurred in a very long time. These 
galaxies may harbor a thick disk in their outer region, while 
the thin disk was stripped early. However, until the struc- 
tural components that make up these outer regions are re- 
solved (e.g. disk versus spheroidal), we can only speculate 
about the nature of these galaxies. 



9 CONCLUSION 

We have presented optical (SDSS g and r) and NIR (J 
and/or H) surface photometry for a sample of 59 SO galax- 
ies covering a range in stellar mass and light concentra- 
tion. Radial age and metallicity gradients out to at least 
5 R e are derived from comparison of the observed g—r and 
r— H(and/orr-J) colors with SP models. 

We find an average central light-weighted age of ~ 4 Gyr 
and central metallicity of [Z/H] ~ 0.2 dex. For most of the 
galaxies in our sample, a large negative metallicity gradient 
is found with radius, with an average of A[Z/H]/Alog(r/R e ) 
= -0.6 ± 0.5. Radial age trends are more heterogeneous. An 
increase in age with radius is found for 58% of our sample, 
a decrease for 19%, and little change for 23% of it. 

The outer regions of 14 out of 59 SO galaxies have very 
old light- weighted ages (> 10 Gyr) and also exhibit large in- 
creases in light- weighted age from the center outward. These 
old outskirts galaxies are found in a range of environments, 
masses, and concentrations. None of these galaxies shows 
nascent spiral structure. Determining the structural compo- 
nent (i.e. disk versus spheroidal) that makes up the outer 
region is a necessary step to further probe the formation 
history of these galaxies. Assuming that the outer region is 
a disk, old outskirts galaxies may be showing the thick disk 
in the outer region, while the thin disk was stripped early. 
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It is possible that old outskirts galaxies require no major 
merger to have occurred in a very long time. 

We find that mean light-weighted age and metallicity 
correlate with both mass and concentration; galaxies with 
both lower mass and lower concentration have, on aver- 
age, younger central ages and lower metallicities than other 
systems in our sample. The mass-metallicity relation ex- 
tends to the outer regions of the galaxy. The observed mass- 
metallicity trends in SO galaxies are consistent with a merg- 
ing formation of these galaxies. The combination of a sta- 
tistically incomplete sample regarding environment (e.g., all 
of the observed lowest mass galaxies are in the Virgo clus- 
ter high density environment) and a scarcity of high density 
environment galaxies makes it difficult to draw conclusions 
regarding the effect of environment on age or metallicity. 
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Table 1. Galaxy Characteristics 



Galaxy 


Morph* 


a? 


bi 


v r i 


Distance^ 


Inch* 


Type 






(') 


(') 


( km s — 1 ) 


(Mpc) 


(°) 




UGC04330 


SBO 


1.1 


1.0 


4865 


69.4 


34.6 


B 


UGC04596 


SO 


1.2 


1.0 


9439 


131.7 


25.4 


T 


UGC04599 


SO 


2.1 


2.1 


2072 


31.8 


18.3 


T 


UGC04631 


so 


1.1 


1.0 


4159 


60.6 


24.3 


N 


UGC04639 


so? 


1.4 


1.3 


8556 


120.0 


23.3 


N 


UGC04737 


so? 


0.8 


0.5 


3813 


56.0 


53.3 


N 


UGC04869 


so? 


2.0 


0.7 


6889 


97.8 


62.9 


T,B 


UGC04901 


SO-A? 


1.1 


1.1 


8424 


118.7 


28.9 


N 


UGC04910 


so 


1.0 


0.6 


8353 


117.8 


28.9 


N 


UGC04916 


so 


1.2 


0.9 


8785 


123.6 


36.2 


N 



The full table appears in the electronic edition of the Journal. 



* From the Uppsala General Catalog (UGC) unless otherwise noted 

1 Major and minor diameter, and radial velocity, from the NASA/IPEC Extragalactic Database (NED) 
k Distances are corrected for Virgo flow and the Great Attractor, from NED 
' T (transition galaxy), B (barred, from Hyperleda), N (normal) 
m From NED if morphology; from Hyperleda if inclination 



Table 2. UH 2.2-m Observations 



Galaxy 


T(sec^ - 


date 


Z.P.(mag)^ 




NS,£ 


UGC04330J 


480 


Apr 2007 


-1.916 


0.037 


10 


UGC04596H 


480 


Mar 2008 


-1.731 


0.004 


3 


UGC04599J 


480 


Apr 2007 


-1.713 


0.034 


11 


UGC04631J 


480 


Apr 2007 


-1.707 


0.037 


7 


UGC04639J 


480 


Apr 2007 


-1.628 


0.030 


9 


UGC04737H 


480 


Mar 2008 


-1.819 


0.201 


10 


UGC04869H 


480 


Mar 2008 


-1.853 


0.056 





UGC04901J 


480 


Apr 2007 


-1.596 


0.025 


9 


UGC04910J 


480 


Apr 2007 


-1.579 


0.057 


7 


UGC04916J 


480 


Apr 2007 


-1.694 


0.018 


7 



The full table appears in the electronic edition of the Journal. 



a Total exposure time in seconds 
^ Zero-point magnitude 

c Standard deviation in the derived zero-points for that field 
^ Number of standard stars in the field 



Table 3. Derived Radial Quantities 



Inner Outer 



Name 


cutoff^ (") 


age(Gyr) 




-tj 


[Z/H] 




-a 


age(Gyr) 


+cr 




[Z/H] 


+cr 


-<7 


H band derived 




























UGC04596 


38 


10.53 


0.81 


0.32 


0.06 


0.01 


0.00 


2.49 


0.15 


0.14 


-0.47 


0.03 


0.04 


UGC04737 


62 


4.69 


0.20 


0.18 


0.22 


0.01 


0.01 


7.31 


1.57 


1.13 


-0.22 


0.05 


0.06 


UGC04869 


100 


3.97 






0.50 






2.88 


0.72 


0.25 


0.03 


0.05 


0.05 


UGC05094 


60 


3.48 


0.15 




0.48 


0.02 




5.51 


3.64 


2.19 


0.04 


0.07 


0.06 


UGC05182 


50 


5.83 


0.43 


0.21 


0.23 


0.01 


0.01 


13.19 




2.35 


-0.50 




0.05 


J band derived 




























UGC04330 


43 


2.08 






0.50 






3.00 


0.77 


0.21 


0.07 


0.05 


0.05 


UGC04599 


58 


2.73 


0.05 


0.06 


-0.01 


0.02 


0.02 


1.05 


0.02 


0.02 


-0.48 


0.02 


0.02 


UGC04631 


3(5 


3.90 


0.42 


0.16 


0.19 


0.02 


0.01 


13.69 




1.08 


-0.69 




0.23 


UGC04639 


58 


4.76 


0.27 


0.24 


0.28 


0.02 


0.02 


2.17 


0.33 


0.25 


0.19 


0.03 


0.04 


UGC04901 


90 


3.49 


0.08 


0.03 


0.40 


0.01 


0.01 


2.50 






0.47 






'he full table appears in the electronic edition of the 


Journal. 



















a Outer radial cutoff in arcseconds 



Table 4. Derived Total Quantities 



Name 


rtot 


gtot 


Htot Jtot 


Re r (") 


g-rtot 


r-H 


UGC04596 


14.0 


14.7 


11.4 


9.9 


0.7 


2.6 


UGC04737 


13.2 


14.0 


10.4 


8.2 


0.8 


2.8 


UGC04869 


12.9 


13.8 


10.0 


11.9 


0.8 


3.0 


UGC05094 


14.0 


14.8 


11.2 


8.3 


0.8 


2.8 


UGC05182 


13.0 


13.8 


10.4 


10.6 


0.9 


2.6 


UGC05403 


13.3 


14.1 


10.6 


10.7 


0.8 


2.7 


UGC05419 


12.3 


13.0 


9.7 


19.0 


0.7 


2.6 


UGC05568 


12.1 


12.8 


9.7 


15.3 


0.7 


2.4 


UGC05766 


12.2 


12.9 


9.6 


17.2 


0.8 


2.6 


UGC08886 


12.9 


13.7 


10.1 


6.1 


0.8 


2.8 



The full table appears in the electronic edition of the Journal. 



Table 5. Binning Scheme 



Bin 


Radial range 


1 


0.0 <r/R e < 0.5 


2 


0.5 <r/R e < 1.5 


3 


1.5 <r/R e < 2.5 


4 


2.5 <r/R e < 3.5 


5 


3.5 <r/R e < 4.5 


6 


4.5 <r/R e < 5.5 



Table 6. Derived Global Properties 



Name 


Concentration 


Mass 


Local Density 






(M ) 


(log Mpc~ 3 ) 


UGC04330 


4.4 


1.22c+ll 


0.40 


UGC04596 


5.3 


1.02c+ll 


-1.41 


UGC04599 


4.4 


6.97c+09 


0.29 


UGC04631 


4.9 


3.88c+10 


0.45 


UGC04639 


4.7 


2.29c+ll 


-1.22 


UGC04737 


4.7 


4.92c+10 


-0.07 


UGC04869 


4.7 


1.87c+ll 


-0.29 


UGC04901 


4.2 


3.62c+ll 


0.43 


UGC04910 


4.9 


3.24c+ll 


0.94 


UGC04916 


5.3 


1.95e+ll 


-0.72 



The full table appears in the electronic edition of the Journal. 



Table 7. Results from Statistical Analysis 



Groups Result 



Age Trends 



inner region ages from J-band and H-band for galaxies in common 





81 


outer region ages from J-band and H-band for galaxies in common 





81 


inner ages of high and low concentration galaxies 





03 


outer ages of high and low concentration galaxies 





20 


inner ages of high and low mass galaxies 





03 


outer ages of high and low mass galaxies 





06 


ages for all bins of high and low concentration galaxies 


1 


90E- 


ages for all bins of high and low mass galaxies 


1 


70E- 


ages for all bins of high and low mass galaxies with low concentration 





09 


ages for all bins of high and low mass galaxies with high concentration 





02 


ages for all bins of high and low concentration galaxies with low mass 





19 


ages for all bins of high and low concentration galaxies with high mass 





04 


age differences of transition and featureless disk galaxies 


1 


20E- 


outer ages of high and low concentration non-old-outskirts galaxies 





34 


outer ages of high and low mass non-old-outskirts galaxies 





05 


outer ages of high and low density non-old-outskirts galaxies 





99 


mass of old-outskirts and non-old-outskirts galaxies 





71 


concentration of old-outskirts and non-old-outskirts galaxies 





91 


environmental density of old-outskirts and non-old-outskirts galaxies 





36 


Metallicity Trends 






inner region metallicity from J-band and H-band for galaxies in common 





81 


outer region metallicity from J-band and H-band for galaxies in common 





81 


inner metallicity of high and low concentration galaxies 


7 


49E- 


outer metallicity of high and low concentration galaxies 





10 


inner metallicity of high and low mass galaxies 


6 


77E- 


outer metallicity of high and low mass galaxies 





01 


metallicity for all bins of high and low concentration galaxies 





09 


metallicity for all bins of high and low mass galaxies 


8 


10E- 


metallicity for all bins of high and low mass galaxies with low concentration 


5 


90E- 


metallicity for all bins of high and low mass galaxies with high concentration 





53 


metallicity for all bins of high and low concentration galaxies with low mass 


9 


1E-4 


metallicity for all bins of high and low concentration galaxies with high mass 





06 


inner metallicity of high and low mass galaxies with low concentration 


2 


5E-4 


inner metallicity of high and low mass galaxies with high concentration 





20 


inner metallicity of high and low concentration galaxies with low mass 


2 


5E-4 


inner metallicity of high and low concentration galaxies with high mass 





68 


metallicity gradient (A[M / H] / Alog(r)) of high and low mass galaxies 





30 


metallicity gradient (A[M / H] / Alog(r)) of high and low concentration galaxies 





30 



Probability from Kolmogorov-Smirnov test that the two groups arc drawn from the same population. Sep- 
aration of concentration groups is at C28 = 4.7, separation of mass is at lxlO 11 Mq, and separation of 
density is at -0.3. 
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Figure 1. H-band radial SB profile for UGC4737. The SB profile 
is shown as black dots with ±<r sky error envelopes as red lines. 
The sky error envelopes are calculated using sky values adjusted 
to sky = sky or i g ±la ays . Black error bars are statistical SB errors 
from profile fitting. The upper axis indicates the radial extent 
from the center in arcseconds while the lower axis is scaled by the 
r-band half-light radius. The outer cutoff is shown as a vertical 
blue dashed line. 
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Figure 2. Radial SB profiles in H-band (black), r-band (ma- 
genta), and p-band (blue) for UGC 4869. The vertical dashed 
cyan lines indicate the separation between the inner and outer 
regions. 1 a statistical SB errors are shown at each radial point. 



Figure 3. Radial SB profiles in g-band (blue), r- 
band (magenta), J-band (red), and H-band (black) 
are shown. UGC or VCC numbers of the galaxies 
are given in the upper right corner of each panel. 
[See the electronic edition of the Journal for the remainder of our sample.] 
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Figure 5. Left: Histogram of total stellar masses for 59 galaxies. 
Middle: Histogram of the galaxy light concentration, C281 f° r 59 
galaxies. Right: Histogram of environmental number densities (in 
log Mpc — 3 ) for 59 galaxies. 



Figure 4. Color profiles in r-H (red) and g-r (blue) are 
shown. The lower axis indicates the radius in terms of the 
r-band half-light radius (R e ) while the upper axis is in units of 
arcseconds. The ±1 a error bars are computed in quadrature 
from each of the combined SB profiles. UGC or VCC galaxy 
numbers are shown in the upper right corner of each panel. 
[See the electronic edition of the Journal for the remainder of our sample.] 
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• x 



log(M,/M Q ) 

Figure 6. Total stellar mass versus light concentration (C28) for 
the sample of 59 galaxies. The point style designates local envi- 
ronmental density in Mpc~ 3 (green circles: d > 0.0, red triangles: 
-1.0 < d < 0.0, and black crosses: d < -1.0). 
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Figure 8. r-H vs. g-r color-color diagram for UGC 10391. 
Galaxy colors are overlaid on a CB10 SSP model grid. Model 
lines of constant age (shown as dashed red lines) increase, left to 
right, from 0.8, 1.0, 2.5, 5.0, and 13.8 Gyr. Model lines of constant 
metallicity (shown as dotted blue lines) increase from bottom to 
top from [Z/H] = -2.2, -1.6, -0.6, -0.3, 0.1, and 0.5. Small, filled 
circles are the average colors of the galaxy's radial bins (the bin- 
ning scheme is noted in Table 5) and each bin is connected by 
the solid line. The central binning region is designated by a green 
star. The error bars on each radial bin represent ± 1 <r, where a 
is the standard error in the mean based on the scatter in color 
for each radial point in the designated bin added in quadrature 
to the sky error as discussed in section £|4, 31 A foreground screen 
dust model color vector with A„ = 0.3 is shown in the upper left 
corner. 
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Figure 7. Total color (u - r) versus t he stellar mass, in units o f 
log solar mass. The dotted line from Kannappan ct al. I j2009aT ) 
separates the red and blue sequences of galaxies. Our full sample 
can be thought of as a collection of red sequence galaxies. 




Figure 9. r-H vs. g-r color-color diagram with CB10 models 
for various star formation histories. Overlaid are models using a 
constant star formation history (cyan), an exponentially declin- 
ing star formation history (red), and a simple stellar population 
(black). Both the constant and exponential models assume a max- 
imum age of 13 Gyr. Lines of constant age are dashed and lines of 
constant metallicity are solid. Metallicity increases from bottom 
to top for all 3 models as [Z/H] = -2.2,-1.6,-0.6,-0.3,0.1,0.5. 
Age lines for the SSP model are as in Fig. [8] Age lines for the 
exponential model increase from left to right as time constant 
t = 100, 13,6.5,4.0,3.0, and 0.1. Age lines for the constant star 
formation model go from right to left as time since star formation 
has occurred = 0.1, 10, 12, 12.8, 12.9, and 13 Gyr. 




0.2 0.4 0.6 0.8 



g-r 




0.2 0.4 0.6 O.fl 



g-r 

Figure 10. Color-color diagrams using the BC03 model (in black) 
overlaid on other population synthesis models (in color). Lines 
of constant age and metallicity are represented by dashed and 
solid lines, respectively, and var y in color accord ing to their val- 
ues for all but the BC03 model. nVlarastor] l|2005h and CB10 SSP 
model grids are s hown in the top and bottom panels, respectively. 
iMarastonl lj2005T ) ages inc rease from left to right as 0.8, 2.0, 5.0, 
10.0, and 14.0 Gyr SSP. IMarastonl ||2005| ) metallicity increases 
from bottom to top as [Z/H] = -2.25, -1.35, -0.33, 0.35, and 0.67 
dex. The BC03 and CB10 SSP model grids have ages and metal- 
licities as in Fig. \E\ 
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H band-derived age 

Figure 11. Comparison of ages derived from J- and H- bands for 
the inner (blue circles) and outer radial regions (red crosses) for 
the 11 galaxies observed in both bands. No error bar is shown 
when there is a lack of a measurable error. 
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Figure 12. Color-color diagram for the entire sample of 41 SO 
galaxies with H-band imaging. The model grid is the same as in 
Fig. [8] The error bar in the upper left corner is the same as in 
Fig. \E\ but averaged for each galaxy at each radial bin. 
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Figure 13. Distribution of average age (left) and metallicity 
(right) of the central 0.5 R e for the full H-band sample, top, 
and J-band sample, bottom. Red lines show only galaxies that 
are observed in both the J- and H-bands. 



Figure 15. Top: Distribution of light-weighted age differences (in 
Gyr) from the inner to outer regions (outer age - inner age) for all 
galaxies in our sample. Middle: Distribution of age differences for 
featureless disk galaxies. Bottom: Distribution of age differences 
for transition galaxies. 




5 10 
Age (Gyr) 



Figure 14. Top and bottom panels show the distribution of av- 
erage ages (in Gyr) of the inner (R < 0.8 Re) and outer (R > 1.2 
R e ) regions, respectively, for all galaxies in our sample. 
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Figure 16. Color-color diagrams for the sample of old outskirts 
galaxies in the H-band (top three rows) and J-band (bottom 3 
rows) plot individually in a panel. The model grid and errors are 
the same as in Fig. [8] 




0.2 



0.4 0.6 

fi — r 



2 - 




0.2 



OA 



0.6 



0.P, 



Figure 17. Color-color diagrams for the low (red) and high (blue) 
stellar mass groups. The same model grid is used as in Fig. [8] the 
top panel shows r-H vs g-r for the H-band sample and the bottom 
panel shows r-J vs g-r for the J-band sample. The red and blue 
solid lines delineate galaxies with mass < 1 X 10 11 Mq and mass 
> 1 X 10 11 Mq, respectively. The solid colored error bars on each 
radial bin represent the standard error in the mean based on the 
scatter in color among galaxies combined in quadrature with the 
standard error in the mean of the systematic sky error. The black 
dotted error bar on the first radial bin represents the galaxy-to- 
galaxy scatter. The reddening line, ages and metallicities, and 
symbols are the same as in Fig. \E\ 




Figure 18. Color-color diagrams for galaxies separated by high Figure 19. Age versus stellar mass (top panel) and concentration 

and low light concentration. The solid red and blue lines represent (bottom panel) for inner (solid circle) and outer (open triangle) 

the average colors for galaxies with C28 < 4.7 and C28 > 4.7, radial regions. Colors denote the radial trend; red indicates an 

respectively. The model grid and error bars are as in Fig. 1171 increase in age from the inner to the outer region, blue indicates 

a decrease in age, and black indicates AAge < 1 Gyr. Arrows on 
the error bars indicate the lack of a measurable error. 
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Figure 20. Color-color diagrams for high (top) and low (bottom) 
light concentration galaxies, separated by C28 = 4.7 are shown 
for the H-band sample. The solid red and blue line represents the 
average colors for galaxies with mass < 1 X 10 11 Mq and mass 
> 1 X 10 11 Mq , respectively. The model grid and error bars are 
as in Fig. [17] 
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Figure 21. Color-color diagrams for high (top) and low (bottom) 
mass galaxies, separated by a mass of 1 X 10 11 Mq, are shown 
for the H-band sample. The solid red and blue line represents 
the average colors for galaxies with C28 < 4.7 and C28 > 4.7, 
respectively. The model grid and error bars are as in Fig. 1171 





Figure 22. Metallicity versus stellar mass (top panel) and con- 
centration (bottom panel) for inner (solid circle) and outer (open 
triangle) radial regions. Colors denote the radial trend; blue indi- 
cates a decrease in metallicity from the inner to the outer region 
and black indicates A [Z/H] < 0.25. Arrows on the error bars 
indicate the lack of a measurable error. 
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